Mixing and the Jet Flap 


B. S. STRATFORD 
(National Gas Turbine Establishment) 


SUMMARY: One of the main problems associated with the “ Jet Flap ” concerns 
the discrepancy in thrust between idealised theory and the experimental results. 
This discrepancy is attributable to the mixing with the surrounding flow of the 
thin two-dimensional jet while still in close proximity to the aerofoil. The 
effect of the mixing may be calculated to a first approximation from a formula 
derivable from first principles, while certain second order effects, which can 
be significant, may be considered qualitatively. 
It is concluded that 
(i) the full thrust* should be experienced by a jet flapped aircraft at cruise, 
(ii) it should be possible to attain a low form drag at cruise in comparison 
with a conventional aircraft, 


(iii) at take-off, an aircraft using a shrouded jet flap would have better thrust 
recovery than one using a pure jet flap (which shows appreciable losses), 


(iv) the use of by-pass engines would further improve the thrust characteristics, 


(v) the practical gains from thrust augmentation, as obtained by controlling 
the mixing, seem likely to be small. 


1. Introduction 


In the historical sketch”? which formed the first paper in this series on the 
“Jet Flap ” an outline was given of the mechanism whereby lift was induced. The 
“Thrust Hypothesis ” (according to which the forward thrust on a jet flapped aero- 
foil remains constant independently of increase in the deflection angle of the jet) also 
was discussed, but only incompletely, consideration being limited to an inviscid flow 
having no mixing. Extension of this hypothesis to a real flow is of interest not only 
because of the scientific curiosity attached to the hypothesis but also because, as has 
been indicated in Part II of the recent lecture®, its realisation might be of some 
economic importance. 


*The “full ti:rust” is here defined as the thrust which would be experienced by that conventional 

aircraft for which the effects of the mixing of the jet-stream may be assumed negligible and 
for which certain conditions at the nozzle, namely mass flow, total pressure, and total 
temperature, are the same as for the jet flapped aircraft. 
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Variation of aerofoil thrust with jet deflection angle. 


Of the experimental evidence so far available, typical values for the ratio 
(forward thrust)/(total jet reaction), taken from Dimmock’s measurements®: “”, are 
shown in Table I and Fig. 1. (These values include losses due to secondary flow 
from the end wall boundary layers and, for the 90° test, due to separation of the 
model boundary layer, but the trend indicated is probably correct). It will be seen 
that, for an angle of deflection greater than about 30°, the forward thrust 
experienced by the aerofoil is greater than the direct forward component of the 
reaction; there are, however, for all values of 6, appreciable losses as compared with 
the prediction of the hypothesis. 


Ramifications of the experimental results are apparent. For an aircraft at 
cruise, when probably the same jet system as at low speeds would be in operation. 
the jet deflection would be sensibly zero and the results of Table I would suggest a 
loss of between 5 and 10 per cent. of the cruising thrust—sufficient to influence 
seriously the economy of the aircraft. For take-off, and to some extent at landing, 
the effects of a loss of thrust would include a reduction in the usable lift coefficient. 
This is because the maximum lift that may be used on an aircraft equipped with the 
jet flap is most probably determined by the thrust that is available for overcoming 
the induced drag. 
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MIXING AND THE JET FLAP 


TABLE I 


VARIATION OF AEROFOIL FORWARD THRUST WITH JET DEFLECTION ANGLE 


Jet deflection angle 6 0° 90° 
forward thrust 

0-85 0-65 0-30 
total jet reaction 

0:529 0 


cos @ 1:00 0°855 


Notes: (1) The values have been corrected to allow for the direct skin-friction drag. 
(2) The results quoted typify ranges of C, corresponding to cruising conditions for 
the 0° test, and to take-off for the remainder of the tests. 


The problems concerning thrust may therefore be put as follows : — 

(i) What is the cause of the loss of thrust as shown in Table I? 

(ii) Can the thrust loss, if any, be predicted for aircraft cruising conditions? 
(iii) Can the loss be predicted for take-off and landing? 


There is also the further problem 


(iv) Might it be feasible so to combine the jet-aerofoil-flow system as to 
obtain more thrust than is predicted by the hypothesis? 


In answer to the first problem it is suggested that the cause of the thrust loss is 
the mixing which occurs between the jet and the main stream. This suggestion 
forms the basis of the present paper and it is found convenient to use three methods 
for its investigation. The main derivation uses energy considerations, the argument 
being that the useful power expended on the aircraft for a given engine power must, 
in some manner, be changed by an amount that will just counterbalance any change 
in the power dissipated by the mixing. The analysis leads to a simple algebraic 
formula and one of the points demonstrated is that the density of the jet has an 
important influence on the thrust in cruising flight. The second and third methods 
for investigating the mixing respectively use momentum principles and consideration 
of the sink effect of the jet in absorbing the external stream. The sink effect 
provides a simple physical picture for the flow, showing, for example, how the 
increment of thrust or drag is induced on the aerofoil. A later paper will include 
both the longer analysis based on momentum, as this is perhaps more readily 
acceptable than are the energy considerations, and also some further details 
concerning the sink effect. 


The theory secures broad agreement with the experimental facts: nevertheless 
in several types of flow it is still unable to provide a satisfactory prediction. This, 
it is suggested, is because the surrounding air has a velocity component 
perpendicular to the jet (in order to reach and mix with it) whereas the basic 
analysis assumes that the jet and main stream velocities are locally parallel. Thus 
there is so far only a one-dimensional theory for the flow locally in the mixing 
region and the effect predicted is accordingly denoted the “ primary” change of 
thrust. An analysis of the second order effects is highly desirable. 
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When the understanding is sufficiently advanced, aerofoil development might 
perhaps be assisted by a semi-empirical approach, such as that used by Squire and 
Young for their analysis of the mixing in a wake. 


net thrust on a conventional aerofoil 

net thrust on a jet flapped aerofoil 

static pressure in the undisturbed main stream 
static pressure at the position of the mixing 
velocity of the undisturbed main stream 


main stream velocity corresponding to p, 


= U;, jet velocity at pressure p, before mixing 


jet velocity at pressure p, before mixing 

velocity difference between the jet and the local main stream 
density, or densities, of main stream fluid 

density of jet fluid 


mass flow rate of the jet 


velocity component of the external flow perpendicular to the 
boundary of the jet mixing region 


ratio of the mass flow rates of the main stream and of the jet in 
the mixing process at pressure p, 


distance measured perpendicular to the jet or boundary layer 
local velocity in the jet or boundary layer 


| (1 - ii) dy, displacement thickness of the bounday layer or jet 


0 
aerofoil chord 


velocity difference between the jet and the main stream at 
pressure p, 


= Ua.,, Velocity difference between the jet and the main stream at 


pressure p, 


refers to conditions at the position of mixing 


refers to conditions in the undisturbed main stream 


refers to the jet 
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and 
In conventional jet propulsion the dimensions of the jet are such that most of 


the mixing occurs so far downstream of the wing as not to affect its pressure 
distribution; hence the thrust is just that of the engine system. The pressure dis- 
tribution might, however, be affected on a jet flapped aircraft, since the mixing of 
its thin jet will occur very close to the wing. This could directly influence the total 
thrust, a possibility which seems the more likely on recalling that when the jet is 
deflected a part of the thrust must be transmitted by the external pressure forces. 
Before looking for a detailed mechanism whereby the thrust might be altered, an 
examination will be made of the energy balance. 


Of the power output from an aircraft jet engine, a part is expended on the 
aircraft and the remainder is eventually dissipated in the mixing between the jet 
and the surrounding air. Thus the energy balance may be written 


Useful power expended on aircraft = Engine power — power dissipated in jet mixing* 


(1) 


Since the useful power expended on the aircraft is equal to the product of its 
thrust and its forward speed, equation (1) gives 


Engine power — power dissipated in jet mixing 
Thrust on aircraft = (2) 
aircraft speed 


he 


in 
The conventional aircraft may be taken as a standard and the jet flap system of 
equal engine power compared with it. Equation (2) then shows that 
Cheech tit — (change of power dissipated in jet mixing) (3) 
aircraft speed 
Thus, if the thrust acting on the conventional aircraft is T., that on the jet flapped 
aircraft 7;, and if the aircraft speed is U,, the proportional change of thrust is 
_ (change of power dissipated in jet mixing) (4) 
T. T.U, 
at 


It is now necessary to examine the mixing process in more detail. 


*Equation (1) assumes that the thermodynamic reconversion to useful work of the power 
dissipated is sufficiently small to be neglected. It also assumes for simplicity that both the 
wake from the aircraft boundary layers and the three-dimensional downwash remain separate 
from the jet. The second assumption does not affect the principle of the argument, since 
the force to be overcome could be external to the aerodynamic system, e.g. gravity during a 
climb. The first assumption appears unimportant provided that none of the speeds are 
appreciably supersonic, 
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For the jet flapped aircraft the initial stages of the mixing take place within the 
pressure field of the wing, a process which, in its essence, may be represented by 
uniform mixing of the jet with 1 times its own mass flow in a region where the static 
pressure is p,, while the remainder of the mixing occurs at the undisturbed pressure 
p,. To begin with, this will be simplified by supposing that it is the whole of the 
mixing which occurs at the pressure p,; then the change in the power dissipation 
from that of the conventional aircraft (where sensibly all the mixing occurs down- 
stream at'the pressure p,) is that due to the change (p, - p,) in the static pressure of 
the mixing region. It will be assumed initially that (p, - p,) is sufficiently small for 
the use of differentials 5p and 6u. It will also be assumed that the velocities of the 
jet and main stream are parallel locally in the mixing region. 


Let the velocity difference between the jet and the main stream be uw. If the 
mixing is complete, that is if the jet mixes with an infinite mass flow, the whole of 
the jet energy based on this velocity difference u, will be dissipated.* Thus, for a 
mass flow m in the jet, 


power dissipation = 4mu,? (whenn—>oo). 


If the velocity difference u, changes by an amount 6u,, due to the change in the static 
pressure, the change in the power dissipation is given by differentiating 
equation (5):— 


change in power dissipation = muta (6u4) (when n —> 00). ‘ . (6) 


Since the change 61, is supposed small, the value of ua in equation (6) may be taken 
either as that before or that after the change. By choosing the value for the 
conventional aircraft, i.e. ua,, (=U; -U,), the velocity difference becomes equal to 
the velocity increase of the jet in passing through the engine system. Thus mu, in 
equation (6) may be put equal to the thrust : — 


Mug & Mua, ‘ ‘ . (7) 
and then, 
change in power dissipation = 7. (814) ‘ (8) 
Equation (4) therefore reduces to 


This relation, whereby the thrust increases if the velocity difference decreases, (and 
vice versa), is fundamental to the present discussion, so that the quantity dy. 
becomes one of the key considerations in estimating the performance. 


*This may be taken as intuitively obvious, by considering axes fixed relative to the main 
stream, or it may be proved from first principles—for example by means of the result that. 
when two masses m, and m, having velocities u , and u, coalesce, the loss of kinetic energy is 


{4m, (u, u,)?}/{1 +(m,/m,)}. 
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If the more complete picture is considered whereby the jet mixes uniformly 


| with n times its own mass flow at pressure p, and with the remainder of the main 


stream at pressure p,, the conclusion is that 


dua (10) 


where U; is the inviscid jet velocity relative to the aircraft in a position where the 
pressure is p,. The value of d5uq is readily calculated on inviscid flow principles 
as follows. 

Take axes relative to the aircraft; let U, be the velocity of the main stream 
where the pressure is p, and let p,; and p, be the densities of the jet and of the main 
stream respectively. Now for any inviscid flow, whether compressible 
or incompressible, 


This follows directly from momentum principles or, for incompressible flow, by 
differentiating the Bernoulli equation of energy, p+4pu?=constant along 
a streamline. 

bp 

pu’ 
Equation (12) illustrates the second feature which is fundamental to the present 
discussion; for a given change in pressure, the change in velocity is small for a flow 
of high velocity or density and large for one of low velocity or density, the product 
pu being the actual criterion. From equation (12) the values of this product pu for 
the jet and for the main stream will determine their respective changes of velocity, 
and hence the change in their velocity difference, and as already shown, this will 
determine the thrust. 


(12) 


Hence, - 


Thus = 8U; BU main stream 


main stream 


The pressure change 5p is the same for the jet as for the main stream with which 
it mixes and so, using equations (12) and (13), 


Thus, from equation (10), 


(15) 
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B. S. STRATFORD 


This result gives the primary change of thrust when a portion of the mixing 
occurs close to the aerofoil and the remainder downstream, instead of the whole of 
it occurring downstream as effectively happens for a conventional aircraft. 


3. Discussion 
3.1. SOME FEATURES AND LIMITATIONS 


(i) In low speed flight, when U,/U; is sufficiently small for the factor 
[1 - (e,U,)/(pjU;)] to be positive, a gain is achieved by allowing the mixing to occur 
where the pressure is less than the free stream pressure. (The action here may be 
likened to that of an ejector pump or thrust augmentor, the external aircraft surfces 
acting in the same way as the internal pump surfaces). On the other hand, a loss 
would result from mixing which occurs where the pressure is higher than that of the 
free stream. 


(ii) If the velocities and densities are such that 


poU,=p\U; (a possible condition) . . (16) 


the primary thrust or drag increment is zero, whatever the position of the mixing. 


(iii) By multiplication by nU,/U;, it will be seen that, when nU,/U; is small, 
the primary change in thrust is proportional to the amount of the mixing which 
occurs at pressure p,. More specifically it is proportional to the momentum 
ratio nU,/U;. 


(iv) For large amounts of mixing the primary change is asymptotic to 
the value 


(v) Algebraically the formula of equation (15) is only a first approximation 
to the equations representing the flow model that was initially chosen. Moreover, 
physically, the chosen model was only a first approximation to the real flow. Now 
the algebraic approximations are tantamount to the assumption that the quantity 
(U,/U, - 1) is small, and Ref. 5 shows that this does not introduce a serious error. 
In contrast the physical approximation is fundamental: it has been assumed that the 
velocities of the jet and main stream are precisely parallel. In the real flow the 
main stream is drawn into the jet, so that the two fluids are not travelling in exactly 
the same direction—a condition which seems likely to cause a superposed loss! 
This “secondary loss” is not considered quantitatively but will be assumed to 
increase with increase of the entrainment angle. 


(vi) Kinetic heating has been neglected in the statement of equation (1) and in 


the derivation of equation (10), but rough calculation suggests that up to a Mach 
number of unity the effect will be small. 
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MIXING AND THE JET FLAP 


3.2. THE VALUE OF n 


To complete the foregoing analysis for the primary increment of thrust, it is 
necessary to find the rate of mixing and hence the value of the mass flow ratio n. 


The standard theory for the mixing between a jet and “still” air has been 
developed semi-empirically. During the mixing the jet “ entrains ” the external air, 
which must therefore have a velocity perpendicular to the jet boundary. (It is 
effectively this component that is neglected elsewhere in the present treatment). 
Denoting this entrainment velocity by 7, it may be deduced from the standard work 
that v is given approximately by 


where U; is the jet velocity or the peak jet velocity. By integrating 7 along both 
sides of the jet the total flow into the jet, and hence the value of n, can be calculated. 
In equation (18) and the subsequent equations (19) to (21) only, both v and U; 
represent local values, which may be deduced from the values at the beginning of 
the mixing region. 


For a jet in moving air, a theory exists“ but no experimental work is known. 
The entrainment velocity as given by the theory can be shown to be satisfied very 
roughly by the relation 


U; U, 


where U, is the velocity of the external air. Since equation (19) appears a highly 
reasonable generalisation of equation (18), it was decided that a crude experimental 
check would suffice for its confirmation. A comparison of the theoretical and 
experimental boundaries for a jet of which the initial velocity was just under five 
times the external main stream velocity is shown in Fig. 2. In the first part of the 
mixing region there is a discrepancy but, as will be mentioned later, this is probably 
attributable to the high turbulence level in the oncoming boundary layer. Thereafter 
the agreement is sufficiently close to justify the use of equation (19) for present 
requirements. Again, integration of v7 on both sides of the jet leads to a value for n. 


No standard work is known concerning the effects of density, but an experiment 
which will be described presently has provided incidental information. It was 
noticed that when a hydrogen jet issued into otherwise still air the suction pressure 
created in the neighbourhood of the jet, i.e. pv”, was roughly equal to that created 
by a similar jet of cold air of the same dynamic head. Thus equation (18) would 
be replaced by 


U; 
v= 55  (U,=0) 


May 1956 93 


or 

ur 

be 

es 

SS 

he 

g. 

ll, 

h 

0 


B. S. STRATFORD 


EXTERNAL MAIN STREAM FLOW. 


EXPERIMENTAL EDGE OF JET —e— 
[oerTH OF Q INDICATES L J 


POSSIBLE RANGE OF Ye IN. In. 
INTERPRETATION OF EDGE | 
THEORETICAL @OGE OF JET ———— 


Fig. 2. 
A comparison of the experimental and theoretical jet mixing rates. 
Jet nozzle velocity=4-8 
*External boundary layer flow is absorbed by jet in regions AB and A’B’. 


Until more evidence is available a reasonable generalisation of equation (19) would 
therefore be 


When the jet is flowing along a surface the mixing rate is thereby halved and 
correspondingly the entrainment velocity v is integrated along only one side. Any 
reduction in the mixing rate on this one side due to a lower turbulence level in the 
jet—the turbulence must become zero at the wall—is neglected. Also, any initial 
laminar mixing region will be neglected, this being counteracted by the increase in 
mixing caused by the turbulence of the main stream boundary layer. The low 
velocity in the boundary layer may be allowed for simply by replacing U, in 
equation (21) by the appropriate boundary layer velocity at the position 
of entrainment. 


One final factor arises from the assumption in the main theory that the mixing 
of the jet with n times its own mass flow is a uniform process in which the velocity 
profiles are rectangular, both for the jet before mixing and for the jet plus main flow 
after mixing. In fact the profile after mixing is peaked at the centre and has a 
“tail” at each side. Some of the entrained air still remains in the tail and so does 
not significantly contribute to the effect on thrust. Consequently the value for n to 
be used in the thrust calculation should be somewhat less than that obtained from 
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equation (21). A factor of 70 per cent. has been adopted. This was chosen after a 
comparison of certain detailed aspects of the square and peaked jets and it is 
thought to be probably within 10 per cent. of the true value; although intuitively 
reasonable, its derivation is by no means rigorous. 


Typical values of n obtained by integrating equation (21) and applying the 
70 per cent. factor range from 0-5 to 1:5 for cruising conditions and from 0°5 to 2:0 
at take-off. It has been assumed that the nozzle width is equal to one third of one 
per cent. of the chord, and that the mixing occurs on one or both sides of the jet 
over a distance of between 5 and 10 per cent. of the chord. 


3.3. MrxING BETWEEN THE AEROFOIL BOUNDARY LAYERS AND THE JET 


The most significant aspect of the mixing with the aerofoil boundary layers is 
that the form drag of the aerofoil may be removed in cruising flight. Only a 
summary of points is presented as fuller details will be given in the later paper. 


(i) Because of the lower velocities in the boundary layer, the rate of mixing 
will be increased, as given by equation (21). 


(ii) The mixing rate will be further increased by the turbulence of the 
boundary layer. 


(iii) Because of the lower velocities in the boundary layer, the angle at which 
the boundary layer fluid enters the jet will, for a given mixing rate, be 
greater than the angle at which the main stream enters the jet; this will 
probably increase the secondary losses mentioned in Section 3.1. The 
effect would be emphasised by the increases in mixing rate mentioned in 
points (i) and (ii) above. 


(iv) The effect of point (iii) would be further emphasised by a high pressure in 
the mixing region, as the high pressure would reduce the boundary layer 
velocities. On the other hand a suction, i.e. a pressure which is less than 
in the free stream, would increase the boundary layer velocities and 
reduce the secondary losses. 


(v) The previous four points have considered the effects of the boundary 
layer on the rate of entrainment and on the secondary losses. There is 
also an effect on the primary increment of thrust considered in the main 
analysis of Section 2. 


Since the boundary layer velocities are smaller than the main stream 
velocity the boundary layer will react more to changes in pressure than 
would the main stream, as shown by equation (12). The corresponding 
factor for equation (15) is [1 - (p,u)/(p,U;)], which remains positive for 
most of the boundary layer under all conditions of flight. Thus there will 
be an increase of thrust if the boundary layer mixing is allowed to occur 
at a high suction, and this increase would hold for cruise as well 
as at take-off. 
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(vi) The previous effect may be reinterpreted on boundary layer theory as a 
reduction or removal of the form drag : — 


The form drag of an aerofoil, i.e. that part of the two-dimensional drag 
transmitted by the normal pressure forces and not by the skin friction. 
is given by 


1 dp 
Co, form 5* dx dx. 
If the mixing between the jet and the boundary layer occurs at a position 
of high suction, the decrease in displacement thickness that results from 
the mixing will persist for the subsequent pressure rise to the trailing edge. 
Thus the value of the integral and hence of the form drag would be 
reduced. A crude estimate suggests that at cruise this could reduce the 
total drag of the aircraft by at least 5 per cent. 


3.4. VISUALISING THE FLOW 


The reasoning so far has suggested that the mixing can affect the thrust acting 
on an aerofoil, but how it does so has not yet been examined. Again only a 
summary is presented, Ref. 5 considering the mechanisms more completely. 


Consider first a jet issuing from an aerofoil when the main stream fluid is at 
rest other than for the action of the jet. The jet mixes with the surrounding fluid 
and absorbs some of it into itself. To do this the jet must create a suction along its 
boundaries. This creation of a suction pressure may be called the “ sink effect ” 
of the jet mixing. 


If the sink effect occurs aft of the aerofoil it will cause a suction on the rearward 
surfaces and this suction will be greater than that produced on the remoter leading 
surfaces. Consequently there will be a pressure drag, and the net thrust experienced 
by the aerofoil will be less than the thrust from the jet alone. 


When the external fluid is not at rest the sink effect of the mixing changes and 
becomes proportional to [1 -(p,U,)/(p;,,U;,,)] where the suffixes “1” indicate 
values at the local pressure p, before mixing. This follows from a consideration of 
the total stream tube area before and after mixing. Thus, for example, if 
p,U,=p,,,U;,,, it may be shown that the flow occupies the same total cross-sectional 
area before and after mixing; consequently there is zero sink effect, and zero effect 
on the aerofoil thrust. This correlates (to a first approximation) with the 
pPoU,=Pp;,oU;,, condition for zero change in thrust as obtained in the main analysis. 


The thrust increment induced on the aerofoil by a given sink also becomes 
changed when the external flow is not at rest. A sink situated downstream of the 
trailing edge still induces a drag force, but a sink situated on the “ body” of the 
aerofoil induces a thrust force. On the body of the aerofoil the sink tends to 
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reduce the velocities in the trailing edge region and to increase them at the leading 
edge, thus causing a forward force. More strictly the sink exerts a thrust or drag 
on the aerofoil according as the local velocity U, at the position of the sink is 
greater or less than the free stream velocity U,. 


As will be shown in Ref. 5, the foregoing arguments lead to a result which is in 
agreement with equation (15). 


4. Application to the Practical Problems. Low Speed Flight 


4.1. TAKE-OFF AND LANDING 


The experiments®: ” quoted in the Introduction have indicated that a pure jet 
flapped aerofoil suffers an appreciable loss of thrust when the jet is deflected as for 
an aircraft at take-off. As a preliminary to a qualitative examination, equation (15) 
may be simplified to 


since, at take-off, p,U,,/(p,U;) is small, while nU,/U; is not large. Equation (22) 
indicates firstly that there is a gain of thrust from mixing where the local main 
stream velocity U, exceeds the free stream velocity U,, and correspondingly a loss 
where U, is less than U,; and secondly that each gain or loss is proportional to the 
amount of mixing which causes it. 


With the jet deflected downwards the pressure on its lower side becomes greater 
than the undisturbed static pressure p,. and consequently (U,/U,,— 1) is negative, 
while on the upper side the pressure decreases and (U,/U,, - 1) is positive. Thus the 
mixing which occurs on the lower side causes a loss of thrust, while that which 
occurs on the upper causes an increase. Since (U,/U,,— 1) for the upper side will 
exceed [-(U,/U, - 1)] for the lower, a net gain of thrust might have been expected; 
but several factors offset this. The velocity difference between the jet and the main 
stream will be greater on the lower side; thus the amount of mixing there will be 
greater (equation (21)) and the loss will be increased. Moreover the boundary 
layer from the lower aerofoil surface is likely to have just experienced a pressure 
rise, so that its flow might be highly turbulent and its velocity very low. Both factors 
further increase the amount of mixing and therefore the value of n in equation (22), 
and so further increase the loss on the lower side as compared with the gain on the 
upper. Probably more important than these factors, however, is the secondary loss 
mentioned in Section 3.1 but omitted from the main theory. On the lower side 
the entrainment velocity v is comparable with the local main stream velocity U, and 
the experiments show that the main flow may go “ straight into the jet”; near the 
nozzle there can be an angle of 120° between the jet and the flow which it is 
entraining. It is therefore to be expected that the secondary losses would be quite 
severe for a pure jet flapped aerofoil. It is not yet possible, however, to make a 
theoretical prediction of these losses. 
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Two trends seem predictable. Firstly, the use of by-pass engines, which would 
allow lower jet speeds for a given thrust, should decrease the mixing rate and there- 
fore the net primary losses. It should also decrease the entrainment angles and 
hence the secondary losses. Secondly, the shroud of the shrouded jet flap also 
should decrease the. mixing losses. The “shrouded jet flap” is a variant of the 
pure jet flap. The jet issues over a mechanical flap, which may be only a few per 
cent. in chordal extent, and follows around the upper surface of this flap by 
“Coanda” action, thereby being deflected when the mechanical flap is deflected. 
The shrouded jet flap has advantages such as simplicity of jet deflection, increased 
effectiveness at cruise, and reduction of induced drag if an engine fails at take-off or 
landing. Its special problems concern the friction and heat transfer between the jet 
and the flap. The name “ shroud,” however, was coined because of the action of 
the mechanical flap in sheltering the jet from the region of very high pressure in the 
corner, just under the jet exit, where very violent mixing occurs for the pure jet flap. 
It is because of this sheltering action that the losses from the mixing on the under 
side should be reduced. The shroud, or mechanical flap, would also increase the 
region of favourable mixing on the upper side of the jet. Thus, in addition to the 
advantages just quoted, the thrust characteristics at take-off and landing would be 
expected to be better for the shrouded than for the pure jet flapped aerofoil. In 
practice attention would have to be paid to conditions at the leading edge where a 
boundary layer separation could cause significant losses; for all the present 
discussion it is assumed that no boundary layer separation occurs. 


4.2. THRUST AUGMENTATION 


The basic theory of Section 2 has suggested that it would be possible to obtain 
an increase of thrust at take-off by causing the jet to mix with the main stream in a 
region of high suction, i.e. in a region where the pressure is much less than that in 
the free stream and where (U,/U,,- 1) is large and positive. With the shrouded jet 
flap a strong suction exists on the upper surface around the “ hinge ” of the flap so 
that. by positioning the jet nozzle at say 80 to 90 per cent. of the chord, the desired 
arrangement would be obtained. With the mixing on only one side of the jet the 
mass flow ratio n might be expected to reach the value n=1-0 in the “ useful ” 
region, and if the mean value of (U,/U,-1) were 2:0, with U,/U,;=1/15, and 
Py/p;=3-0, the increase of thrust would be given by equation (15) as 


T,-T. _ 2(1-1/5) 


According to a quick calculation the skin friction between the jet and the aerofoil 
would largely counteract any gain of this magnitude. 


Supposing, on the other hand, that very rapid mixing could be achieved in 
the suction region. The limit to the gain is given by equation (15) when the value 
of n becomes infinite, i.e. 


tT. l 1- (when n —> 
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JET FLAP 


MIXING AND THE 


ROW OF SMALL 
“HALF FISH-TAIL 


(0) (0) IN THS REGION 
Fig. 3. 
One method of increasing the rate of mixing; additional losses that may result. 


For the values just given this would give 160 per cent. augmentation; even for, say, 
n=5 the augmentation would be 40 per cent. A mixing rate giving n=5 could be 
achieved by special arrangement of the jet, for example by the scheme illustrated 
in Fig. 3, where a row of “half fish-tail” nozzles is placed along the span. The 
aim in this scheme is to increase the peripheral area between the jet and the main 
stream, and to reduce the damping effect of the wall on the mixing turbulence. Also 
the jet skin friction should be less than for a shrouded jet flap. 


To test the foregoing suggestion, and to find whether other effects might offset 
the gain, an experiment was performed with the model shown in Fig. 4. In the 
experiment it was found that thrust was lost instead of gained over the whole range 
of jet velocities, although the amount of this loss was slightly less than would have 
been predicted for a pure jet flap. 


At a condition representative of take-off a value was achieved for n of 4-0. 
The mean value of (U,/U,- 1) was 1:0 in the mixing region, while U,,/U;=1/20 
and p,/pj=1-0. (The incidence of the aerofoil body was - 4° and the upper surface 
of the flap was at 45° to the aerofoil body; the value of the jet coefficient C; was 1-3). 
With these values, equation (15) would predict the primary increase of thrust to 
be 16 per cent. 


The thrust measured on the balance under the given conditions was 83 per cent. 
of the jet reaction, i.e. a loss of 17 per cent. The loss increased slowly to 20 per 
cent. when C; was altered either way by a factor of about four. With U,, equal to 
zero the measured thrust was 75 per cent. at this incidence. 


In the experiment described, some of the losses could be attributed to skin 
friction and to end effects both from the tunnel wall boundary layers and from the 
clearance between the model and the tunnel walls. The results quoted have not 
been corrected for any of these factors. Nevertheless, flow visualisation revealed 
very little secondary flow in the main stream and it is probably necessary to attribute 
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STATIC HOLES 
IN MIXING REGION 
(5 IN ALL) 


MAIN PROPULSIVE JETS 


~SLOT FOR BOUNDARY LAYER 
CONTROL JET 
(CONTROL JET FOUND UNNECESSARY) 
Fig. 4. 


The model for the thrust augmentation experiment (see also Fig. 3). 
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MIXING AND THE JET FLAP 


SANOWICH LARGE 


SECONDARY yar. 
LOSSES FROM 

+ « JET DRAWN MIXING OF 

ONTO SURFACE JET CACKFLOW 


Fig. 5. 
The loss of thrust when “ sandwich ” air causes the jet to impinge on to a surface. 


a considerable proportion of the losses to secondary effects connected with the 
mixing. Experimental confirmation of this would require either a detailed analysis 
of the jet flow or some arrangement that would eliminate the end effects on the 
model. The mixing downstream of the nominal mixing region would also need 
investigation, although its contribution to the thrust loss seems likely to be small. 


In addition to the secondary loss that has previously been discussed, i.e. that 
due to the angle at which the external flow is drawn into each jet, various other 
sources of loss occur in the type of arrangement examined. The essentially three- 
dimensional flow between neighbouring nozzles may cause a loss. Also there is a 
possibility of secondary flows at the corner at the flap “hinge,” as shown in 
Fig. 3 (b). if the jet total head has not become uniform in a spanwise direction before 
reaching the corner. (For a non-uniform jet the high momentum air would tend to 
ignore the corner, while the low momentum air would curl under in order to replace 
it). An additional possible source of loss occurs where the jet has “ sandwich ” air 
between it and the surface, as shown in Fig. 5. The jet will eventually entrain all 
the sandwich air and impinge on the surface at an appreciable angle, say y. There 
will thus be a direct loss of (1- cosy) times the jet momentum and, due to the 
mixing of the reversed flow associated with this momentum loss, there will be, as 
shown in Fig. 5, subsequent losses that might worsen the situation. In the 
experiment described there would be some sandwich air, owing to the nozzles being 
slightly away from the surface, while any arrangement of a jet blowing across the 
slot of a slotted flap, or of a jet issuing from a nozzle parallel to the surface but 
away from it, perhaps in order to increase the mixing while avoiding maldistribution 
along the span, would be likely to meet this difficulty. 


Thus, while thrust augmentation, i.e. the attainment of an increase in thrust by 
controlling the mixing, at first seems an attractive proposition, there are many 
complicating factors to be taken into account. A worthwhile gain would require 
rapid mixing. It seems likely, though not finally proven, that rapid mixing would 
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inherently produce various types of secondary flow, and these secondary flows 
would reduce the gain. The final thrust from the system should be somewhat 
greater than for a pure jet flapped aerofoil, but unless the details of the flow can 
be improved the final thrust would probably not exceed the 100 per cent. predicted 
by the simple thrust hypothesis. 


5. Application to the Practical Problems. Cruising Flight 
5.1. THE PRIMARY CHANGE OF THRUST 


Under cruising conditions the dominant factor for the primary change of 
thrust is [1 -(,U,)/(»,U;)]. 


It so happens that for a typical current aircraft at cruise the value of p,U; is 


approximately equal to that of p,U, — for example the jet velocity may be double the 
aircraft velocity and the jet density a half of the ambient. Thus the value of 
{1 -(p,U,)/(p;U;)] is close to zero, this being true both for aircraft with pure jet 


engines and for those with by-pass engines. Consequently, whatever the position — 
of the mixing, the primary loss or gain of thrust as indicated by equation (15) also | 


would be close to zero. 


5.2. THE OVERALL RESULT 


In addition to the primary effect, both the secondary loss and also the mixing 
with the boundary layer must be taken into account when estimating the 
overall result. 


The secondary loss should be small at cruise: with the high main stream 
velocity the entrainment velocity v perpendicular to the jet would be low, as from 
equation (21), and consequently the entrainment angle, tan~'(v/U,), would be 
very small. 


As regards the mixing between the jet and the aerofoil boundary layers, the 
discussion of Section 3.3 has suggested that, with care, some gain should be possible. 
This gain would be achieved by allowing the jet to absorb the boundary layer in a 
region of high suction—an arrangement that would be practicable on the shrouded 
jet flap—and the result could be represented as a reduction in the form drag. 


Thus, from the foregoing arguments, the net effect from the three aspects of the 
mixing process should be that an aircraft at cruise will experience effectively the 
full thrust and possibly some gain. 


5.3. THE EFFECT OF JET TEMPERATURE 


The foregoing conclusion is to be contrasted with the results for a pure jet 
flapped aerofoil (at zero deflection) using a jet which has the same temperature and 
density as the main stream. 


When the densities of the jet and of the main stream are equal, the factor 
[1 - (p,U,)/(p,U;)] approaches unity, instead of being zero as in Section 5.1. More- 
over, for the pure jet flap the mixing region is just aft of the trailing edge and so is in 
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Fig. 6. 


Drag/thrust increments under simulated cruising conditions, using a low density (hydrogen) jet. 


a position where U,/U,-1 is negative. There is thus no reason to suppose zero 
loss from the mixing. For example, using U,/U;=}4. U,/U,-1=-9-1, and 
n=1-5, equation (15) indicates a primary loss of 2:2 per cent. Since the hot jet 
would give zero primary loss, it is to be concluded that experiments using a cold air 
jet would be inadequate for simulating cruising conditions for a jet flapped aircraft. 
(It will be noticed that the 2:2 per cent. predicted is still a much smaller loss than the 
6 per cent. quoted from the experiments of Ref. 3, in Table 1. To this 2:2 per cent. 
must be added, however, the secondary loss and the additional effects from the 
boundary layer. The boundary layer effects could cause a significant loss with the 
pure jet flap of the experiment even although on an aircraft a gain should result). 


Since the difference which the theory shows between the behaviour of a hot 
jet and that of a cold one arises from the difference in density, the reasoning could 
be tested by using a cold jet of low density gas—assuming that in fact the tempera- 
ture does not cause other effects. Experiments were therefore performed with a jet 
of hydrogen issuing undeflected from a pure jet flapped aerofoil and the results are 
shown, superficially, in Fig. 6. At face value these provide strong confirmation of 
the theory. The tests were in the nature of a pilot experiment, however, while the 
effects sought were not only very small but also very sensitive to test conditions. 
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These results of Fig. 6 were obtained by pressure plotting, whereas-the results from 
tests using a balance were quite inconclusive, the scatter obscuring any definite 
trend. Thus it would be highly desirable to repeat the experiments using more 
refined techniques. Provisionally, the evidence suggests that the theory is correct 
and that under cruising conditions the mixing should still allow the full thrust 
to be obtained. 


6. Conclusions 


The foregoing examination has demonstrated certain principles governing the 
thrust in a jet flap system, but problems still remain. 


Concerning the loss of thrust at take-off or landing, only general trends may so 
far be suggested. It is expected that the shrouded jet flap would have greater thrust 
recovery than the pure jet flap, and that use of by-pass engines with lower jet 
velocities and smaller mixing rates would further reduce the losses. 


Under cruising conditions the indications are that a jet flapped aircraft should 
experience effectively the full thrust and possibly some gain, the gain coming from 
a reduction in the form drag of the wing. 


It is fairly definite that any significant augmentation of the take-off thrust (by 
controlling the mixing) would require a very high mixing rate. It seems likely, but 
not finally proven, that there are various secondary losses inherent in such rapid 
mixing and that these would counteract much of the augmentation. The attainable 
thrust might well be expected to exceed that for a pure jet flapped aerofoil, but it 
seems unlikely that the prediction of the thrust hypothesis would be surpassed at 
large flap deflection angles unless a radical improvement is possible in the 
augmentation system. 


These conclusions follow from the arguments of the basic theory, with super- 
posed qualitative considerations of secondary losses and boundary layer effects. 
The basic theory shows that the primary change of thrust is given by 


T,-T. _ 
T. [1+U,/(nU,)] 


Whereas this primary change is expected to indicate the general trends correctly, 
the secondary loss, due to the component of velocity of the main flow perpendicular 
and into the jet, may be important and it would be desirable to develop a theory 
for its prediction. 


The manner in which the thrust or drag increments are induced on the 
aerofoil may be visualised by considering the sink effect produced by the jet in its 
absorption of the external flow. 
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Numerical Results for the Longitudinal 
Stability Derivatives and Stability Loci 
of a Low Aspect Ratio Rectangular 

Wing in Supersonic Flow 


LEE WINOGRAD* and JOHN W. MILESt+ 


(University of California) 


1. Introduction 


The longitudinal stability derivatives and stability loci for a single-degree-of- 
freedom, pitching oscillation of a rectangular wing having an effective aspect ratio 
less than unity have been calculated“ from the results of a previous, analytical 
investigation. The results extend those previously available for a rectangular 
wing of effective aspect ratio greater than unity®: ”. 


2. Stability Derivatives 


The stability derivatives Cuas and Cy, are plotted in 
Figs. 1 and 2 as functions of the effective aspect ratio (A =(M? - 1)! (span/chord)) 
for a mid-chord (a= 4) axis of rotation; the Mach number M appears as a separate 
parameter only for Cujo, and Cys. The notation is that of Ref. 2. 


3. Single-Degree-of-Freedom Stability Loci 


The stability loci of single-degree-of-freedom pitching oscillations about an | 


axis x=a (the fraction of total chord aft of the leading edge) may be obtained by 
equating the total damping derivative (Cu, +Cxao/a:) to zero and solving the result- 
ing quadratic equation for a. The pairs of values for a so obtained are plotted as 
a function of M with A as parameter in Fig. 3. The Mach numbers at which the 
roots of this quadratic equation coalesce (i.e., the peaks of the curves of a plotted 
against M) are plotted against true aspect ratio (A’=span/chord) to obtain the 


*Graduate Student, Department of Engineering, Berkeley. 
+Department of Engineering, Los Angeles. 
Received March 1955. 
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Fig. 2. 
Moment coefficient derivatives. 


stability envelope of Fig. 4. These results reveal that single-degree-of-freedom 
instability is impossible for effective aspect ratio less than 0-65 (to two figures), i.e., 


+ Cytaa pat <0 if A<065 . (1) 


4. Accuracy 


The stability locus of Fig. 3 for A=1 is in agreement with the exact (linearised) 
result of Ref. 4 to within 1-2 per cent. This, together with the maximum 
discrepancy of 3 per cent. in centre of pressure at A=1 (cf. Ref. 2), implies that the 
foregoing results, which are based on an expansion in the effective aspect ratio A, 
are adequate for practical application in the range 0 < A <1, 
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Some Integrals Relating to the Vibration of 
a Cantilever Beam and Approximation for 
the Effect of Taper on Overtone Frequencies 


A. I. MARTIN* 


SUMMARY: In the investigation of approximate numerical values of overtone 
frequencies of a turbine blade it is desirable to know formulae for certain 
integrals of the modes of vibration of an ordinary cantilever beam. The first 
object of this paper is to obtain such formulae and to arrange them in 
tabulated form. By proceeding along the lines of the calculus of perturba- 
tions, these results may then be used to obtain new formulae which give 
second order approximations for the effects of uniform breadth and thickness 
tapers on the overtone frequencies. The theory gives good agreement with 
experiments for tapers which do not exceed about 0-5, which is present-day 
practice for turbine blading. 


1. Introduction 


The mathematical solution of the ordinary cantilever beam is well known, the 
frequencies and corresponding modes of vibration being given by elementary 
formulae quoted in the textbooks. In the theoretical investigation of the flexural 
vibrations of a turbine blade, or propeller, e.g. in the investigation of rotary effects 
and twist correction factors, it is often desirable to have, in a definite form, 
formulae for the integrals of these modes. Part I of this paper considers some of 
these integrals. 


If a cantilever beam is uniformly tapered the problem of obtaining its overtone 
frequencies is much more difficult. The solution is exact in a few cases only; 
namely, when the beam has a sharp point or a sharp edge at its free end® * *. 
Several authors have made approximations to the general problem, both by using 
the theory of integral equations and by extending Rayleigh’s device for approxima- 
tion to the fundamental.t Part IT gives second order approximations for the effects 


*The work described here was done while the author was with the Bristol Aeroplane Company. 


tA list of some of the most efficient methods of approximation will be found on pp. 144-152 
of Ref. 1. 
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of uniform breadth and thickness tapers. The formulae may be obtained by 
proceeding along the lines of the calculus of perturbations and by using the integrals 
given in Part I. It is assumed that the frequency of vibration can be expressed as 
a double power series in two taper parameters and that higher powers than two can 
be neglected. Thus the paper does not cover the range of “extremely large ” 
tapers, but concentrates more on the range applicable to present-day turbine blades. 


NOTATION 
y mode of vibration of cantilever beam of length / 
frequency number 
p = tan 47 
I length of the beam 
x distance measured from fixed end 
a,8 breadth and thickness parameters (see Section 5) 


Fy, = V(A/A,), ratio between frequencies of tapered and untapered 
beams 


E —_Young’s modulus 

moment of inertia of cross section 
m= per unit length 
K,,K, see equation (3) 
A,B __ constants defined by equation (4) 
A,,B,,.C,,D, constants defined by equation (24) 
b,t breadth and thickness of beam 
p angular frequency 


X= m,p*/(El,) 


4,.@,,.... Coefficients in the expansion of F,, in powers of z and 8 
k = 2/B 
\ 1s A, \ 
- see equation (18) 
anne see equation (16) 
Vir 


Suffix 0 denotes values when z= 8 =0, i.e. with no taper 
Dashes denote differentiation with respect to x 
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VIBRATION OF A CANTILEVER BEAM 


Part I 


SOME INTEGRALS OF THE MODES OF VIBRATION OF THE ORDINARY CANTILEVER BEAM 


2. Fundamental Relations 
The differential equation for flexural vibration of a cantilever beam is 


with boundary conditions y=y’=0, at x=0, 
and y”=y”’=0, at x=. 


Dashes denote differentiations with respect to x, / is the length of the beam 
which is supposed to be rigidly fixed at x0 and free at x=/, and the dimension- 
less number » is related to the angular frequency p by the formula 


The solution of equations (1) and (2) is 


y=K, (x)- K, (x), (3) 


where K, (x)=A cos +B sin 


K, (x)=A cosh +B sinh 


and A and B are constants to be determined from the pair of equations 


A (cos 7 + cosh n) + B (sin » + sinh n) = 0, (4) 
A (- sin + sinh »)+ B (cos » + cosh n) = 0. 
Elimination of A and B gives the well-known equation 
cos cosh n= - 1; ‘ ‘ (5) 
the first few roots of which are given by Table I. 
TABLE | 
Mode ; Fundamental Ist Overtone ‘2nd Overtone 3rd Overtone 4th Overtone 
of or or or or or 
Vibration Null Harmonic Ist Harmonic 2nd Harmonic 3rd Hermonic 4th Harmonic 
1-875 4-694 7-855 10-996 14137 
p=tan 1:362 —1:019 1-001 — 1:000 1-000 
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The formulae may be greatly simplified by introducing the parameter p= tan 4)). 
giving 

2 2 

and cosh »= (6) 


The last of equations (6), which is obtained from equation (5), shows that p? > 1. 
From the graphs of cos 7 and -1/cosh » it readily follows that p > 1 for the even 
harmonic and p < - 1 for the odd harmonics. Also 


2p 
sinh n='+ 


where the upper sign is taken for the even harmonics and the lower sign for the odd 
harmonics. 


The first of equations (4) now gives B= -(1/p) A for the even harmonics and 
8=pA for the odd harmonics. The constant A, which may be chosen arbitrarily. 
will be taken equal to /-!. 


Then | y*dx=1, as shown later. The formulae for K, (x) and K,. (x) now 
0 
become 


K, (cos sin " *). 


1 


for even harmonics. With odd harmonics p should be replaced by - 1/p. 


In the sections which follow certain integrals involving y and its derivatives 
are given in terms of » and p. The results are presented for the even harmonics: 
for the odd harmonics the corresponding expressions can be derived by replacing 
pby -1/p. 


3. Integrals of K, and K, 
Now K,” (x)= - and K,” (x)= Rm... @ 


Also, by equations (6) and (7), 


K, 0)=K, O)=I"' ?, K,’ (0)= - np } 
K,(O= -K, -I-'"", K,’ (O= - K,’ O= - 
for even harmonics. 


(9) 


Using Green’s formula (uv” — which may easily be 
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verified by integrating by parts, first with u=K,, »=K, and then with u=K,’, 


v= K,’, equations (8) and (9) show that 
l l 
| K, (x) K, (x) dx= K,’ (x) K,’ (x) dx=0 
0 0 


Thus K,, K, and K,’, K.’ form two pairs of orthogonal functions. 


Direct integration of the functions in equation (7) yields 


“ 0 


np 


0 0 


It is also a simple matter to obtain repeated integrals of K, (x)°, etc. 
results are given by 


l F 
4 


= 


and | Ky (dx dx= 1 2). 
"Pp 


In addition, the two formulae 


and 


are required. 
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4. Integrals of y and its Derivatives 


Using equations (3), (8) and (9) it is a simple matter to evaluate y and its 
derivatives at either end of the beam. Thus, in addition to equation (2), 


y °F (14) 


and y= - = atzx=i. 
The integrals given in Table II can be obtained without difficulty by using 


equations (1), (2), (3), (8), (10), (11), (12), (13) and (14). Details of the work are 
given in Appendix I. 


TABLE II* 
1 
(1) | 
0 
(2) yy’dx= 7 
1 
(3) |» dx= ( +2) 
21? 
(4) J» dx= 
1 
Maze 4 2 
(6) | 
0 
1 
4 
(7) y’y"dx= — x 
(8) | 
(9) | y’ 
2 
foram (+ 2) 


*For even harmonics; for odd harmonics replace p by —-1/p. 
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TABLE II (continued) 


3) 


(11) 
0 
l 
(12) | xyy'dx= 4 
2 2 
(13) = =) 
2I\ np 
l 
1 1 
(14) — =) 
l 
1 6 4 
(15) | xo dx = 5) 
F 2 
(16) { xy’y"dx= 
0 
l 
nf 
(17) xy’y""dx= — ap 
2 
2 
t 
(19) xy’ y"dx= — 
2 
(1 4 
(20) | Pdx= (<3 
215 
21) J (3 + + 
2 
l 
3 2 
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By equation (1), the results (1) to (10) in the table determine any integral of the type 


Ax 


in terms of 4 and p, where m and n are positive integers or zero. Similarly, the 
results (11) to (20) yield 


ay ry 
dx" ‘ 


1x. 
The integrals (21), (22) and (23) are required in Part il. 


Part II 


SECOND ORDER APPROXIMATION FOR THE EFFECT OF DOUBLE TAPER ON OVERTON! 
FREQUENCIES 


5. Method of Solution 


By using the integrals in Table II, second order approximations may be 
obtained for the effects of uniform breadth and thickness tapers on the frequencies 
of vibration of a cantilever beam. Let the breadth b and thickness ¢ of the beam 
vary according to the linear relations 


b=b,(1-2%) and t=1,(1-8%), 


where z, 8, b, and f, are constants, / is the length of the beam, and x is measured 
from the fixed end. Thus, the breadth parameter z is equal to the ratio between 
the difference in breadths at the two ends of the beam and the breadth at the fixed 
end; and similarly for the thickness parameter 2. The moment of inertia / and the 
mass per unit length m then vary according to the relations 


Hence, since the ordinary differential equation of flexural vibration is given by 
(Ely”)” =mp*y, where p is the angular frequency. we obtain 


where A=m,p*/(EI,). The boundary conditions are again given by equation (2). 
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It is assumed that A and y can be expanded in powers of z and #3; that is 


+A, Z+A,B+A, 2° +A, 4+... 
and 
where A,, A,, ... are constants and y,, y,.... functions of x. Substituting the values 


of A and y given by (16) in equation (15) and equating like powers of z and £, gives 
a series of differential equations from which the numbers A, A,, ... may be deter- 
mined as simple expressions involving and p. Details of the work are given in 
Appendix II. 


Having thus obtained the numbers A, A,,. . ., the equation F\, = /(A/A,) then 
gives the ratio of the frequency of vibration of the tapered beam to the corres- 
ponding frequency of vibration of the untapered beam. The final result, as far as 
terms of the second order, is given in Table III. 


TABLE III* 


Fy, =1+a,a+a,8 + +... 


_2 
5 1 1 1 
1 19 13 4 4 


Some approximate numerical values are given in Table IV. 


TABLE IV 


Mode Fundamental Ist Overtone 2nd Overtone 3rd Overtone 4th Overtone 
of or or or or or 
Vibration Null Harmonic Ist Harmonic 2nd Harmonic 3rd Harmonic 4th Harmonic 


a, 0°3066 0:0943 0°0324 00163 0°01 
a, 0-114 —0°3115 —0-4353 —0°4673 —0°48 
a, 0°1914 0-051 0-019 0:0099 0:006 
a, 0:0013 —0:0317 - 00052 —0:002 —0-0008 
a 0-077 —0-0419 —0-0593 —0-0548 —0°0575 


*Even harmonics; for odd harmonics replace p by —1/p. 
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6. Examination of the Numerical Results* 


For thickness taper, a comparison has been made between theoretical results 
given in Table III and corresponding experimental values. In the experiment, 
several specimens of various thickness tapers were put into a state of vibration by 
an Electro-Magnetic Induction Exciter, and natural frequencies were obtained by 
resonance: The root thickness of each specimen was } in., the constant breadth 
was 2 in. and the length 6 in. From a knowledge of the material properties, 
density =0-283 Ib./in.* and Young’s modulus=29-2 Ib./in.*, the untapered 
frequencies were calculated. The experimental values of F,, were then obtained 
by dividing the natural frequencies by the corresponding untapered frequencies. 
The theoretical results were obtained by taking F,,=1+a,8+a,6", where a, and a 
have the values given in Table IV. Table V gives the approximate values (z=0). 


TABLE V 


COMPARISON BETWEEN THEORY AND EXPERIMENT FOR THICKNESS TAPER 


Fundamental Ist Overtone 2nd Overtone 3rd Overtone 4th Overton 


Mode of or or or or or 

Vibration Null Harmonic Ist Harmonic 2nd Harmonic 3rd Harmonic 4th Harmonic 

B=0°1 | Theory 1-012 0:9684 0:9559 0:9528 0:9514 

Experiment , 1:0 0:9772 0-9512 0:9450 0:9424 

B=0:2| Theory | 1:026 0:9359 0:9106 0:9044 0:9017 

| Experiment 1:027 0°9487 0:9054 0:9045 0:9298 

B=0°3| Theory 1-041 0-9026 0-8642 0-855 0-8508 

| Experiment 1-054 0-9202 0°8596 0-8515 0°8701 

B=0-4| Theory 1:056 0:8685 0°8165 0°8044 0:7988 

Experiment 1-08 0:9031 0°8179 0-8032 0°8136 

B=0°5 | Theory 1:076 (°8335 0°7677 0°7528 0°7456 

Experiment 1:116 0-8718 0:7752 0°7555 0:7539 

B=06) Theory 1-096 0:7977 0°7177 0-7002 0-6913 

Experiment i132 08063 0:7253 0°7009 0°7005 

B=0°7 Theory 0°7611 06664 06463 06358 

Experiment 1:205 0-7806 0°6775 0°6438 0-6408 


It is also interesting to compare the values given by taking the first six terms 
in the series for Fy, with known theoretically exact values in the case when either 
z or 8 (or both) is equal to unity, i.e. when the cantilever has a sharp edge (or a 
sharp point), at its free end. The results are contained in the tables (a), (b) and (c) 
which follow. 


*The experimental results were kindly supplied by Mr. K. B. Chamberlain of the Bristol 
Aeroplane Company. For the comparison between theory and experiment, the author is 
indebted to Mr. J. Runcie. 
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(a) z=0, B=1; Fx. =1+a,.+a,;. gives the exact values. 


ults Mode Fundamental Ist Overtone 2nd Overtone 3rdOvertone 4th Overtone 
nt, of or or or or or 
by Vibration Null Harmonic Ist Harmonic 2nd Harmonic 3rd Harmonic 4th Harmonic 
by F,,=1+a,+a, 1-191 0°6446 0°5054 0°4779 0°4625 
dth Exact value 
(Ward) 1°512 0°4873 0°4116 
les, 
red 
(b) z=1, B=0; Fy =1+a,+a,. Meyer gives the exact values. 
es. Mode Fundamental Ist Overtone 2nd Overtone 3rdOvertone 4th Overtone 
a of or or or or or 
0) Vibration Null Harmonic Ist Harmonic 2nd Harmonic 3rd Harmonic 4th Harmonic 
F,,=1+a,+a, 1-498 1-145 1051 1-026 1-016 
Exact value 
(Meyer) 2:035 1-409 1-223 1-153 1-118 
(c) «=B=1; Fyr=1l+a,+a,+a,+a,+a,;. Ward’s paper gives the exact 
ne values for a sharply pointed pyramid on a square base. Dr. Dorothy 
Wrinch® has considered a thin conical cantilever bar of circular section. 
nic The two cases agree with one another. 
Mode Fundamental Ist Overtone 2nd Overtone 3rd Overtone 4th Overtone 
of or or or or or 


Vibration Null Harmonic Ist Harmonic 2nd Harmonic 3rd Harmonic 4th Harmonic 


F,,=1+a,+a, 


+a,+a,+a, 1-69 0:7602 0°5516 0°5021 0°4777 

Exact value 

(Ward and 2°48 0:9599 0:6242 ~- 
Wrinch) 


The worst case is (c) in which the first six terms of the series for Fy, may differ by 
nearly as much as 32 per cent. of the total sum. 


As an example in the case of breadth taper, let z=0-5 and 8=0. Equations 
ns (15) and (30) on pp. 146 and 157 of Ref. 1 show that the exact value of F;,, for the 


er fundamental mode of vibration, lies between 1-202 and 1-240. The value obtained 
‘ by taking the first six terms of the series is 1-201, which gives an error of less than 
34 per cent. 


The foregoing considerations show that the results of Table III may be expected 
to give a sufficiently good approximation for tapers which do not exceed 0-5; the 
results appear to be somewhat better for the overtone frequencies than for the 
ol fundamental. As either 2 or 8 approaches unity, the approximation can no longer 
be expected to be satisfactory. 
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Appendix I 
CONSTRUCTION OF TABLE II 


Entries (1) to (10): Integral (1) follows from the first result in equation (10) and the 
first two of equations (11). Integral (2) is immediate by the relation for y in equation 
(14). Integral (3) follows from equations (3), (8) and (11). By integrating by parts, the 
left hand side of integral (4) is equal to 


The equations (3), (10) and the second line of equation (11) yield integral (5). Next, 
integral (6) is immediate by the second line of equation (14), and (7) and (8) follow from 
equations (3), (8) and (11). Integral (9) is immediate by the first equation in the first 
line of equation (14) and, finally, (10) follows from equations (3), (8), (10) and (11). 


Entries (11) to (20): To obtain integral (11), use equation (1) and integrate by parts. 
Thus the left hand side is equal to 


/4 

1) ) 0 


0 0 0 


The first term in the curly brackets is zero and the second is given by integral (5) in the 
table. By equations (3) and (8), the third term is equal to 


t 2 
2 2 2 
x {(K,’)? - (K,’)?} dx= + {(K,’)? - dx - {(K,’)? - (K,’)?} dxdx, 


0 0 


by integrating by parts. Integral (11) now follows from equations (11) and (12). By 
integrating by parts, integral (12) is equal to 


1 1 3 


from equations (2) and (14) and integral (1). 
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Using equations (3) and (8), integral (13) becomes 


l l 
2 

0 0 0 


2 


[? 


(K,? K,”) dxdx 


— 


and the result follows by equations (11) and (12). 


Integrating by parts twice and using equation (2), integral (14) is equal to 


1 
| yy"dx— (P+ 5 (y’)? dx, 


0 0 


and the right hand side is obtained from equation (14) and integrals (3) and (5). For 
integral (15), integrating by parts gives 


l 


[xy] - | | 


and the result follows from equations (2) and (14) and integrals (2) and (13). 


Integral (16) is equal to (I)? - (y’)? dx, 


2 


and the result follows from equation (14) and integral (5). Integral (17) has already 
been considered in the derivation of integral (11). For integral (18), integrating by parts 
gives 


l t 
l 
0 


0 


and these integrals have already been evaluated. 


Integral (19) is equal to minus one half the value of integral (8), by integrating by 
parts. Finally, a further integration by parts shows that integral (20) is equal to 


l 


and the result follows by equation (1) and integrals (9) and (13). 
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Entries (21), (22) and (23): Integration by parts shows that integral (22) is a conse- 
quence of equation (14) and integral (11), and that integral (23) follows from equation 
(14) and integral (21). Lastly, using equation (1), integral (21) is equal to 


l t 
4 1 


l l 
= "dx = { 2 | } 
1) 


0 0 0 


by integrating by parts. The first integral within the brackets is integral (14). By 
equations (3) and (8), the second integral is equal to 


l 
x? ((K,’)? (K,’)"} dx= - [ [cK dx - 


0 0 


0 uv 


by integrating by parts twice. The required result now follows from equations (11), 
(12) and (13). 


Appendix II 


DETERMINATION OF THE COEFFICIENTS A,, A,, A,, Az, A, and A, 


Clearly, on putting = @=0, equation (15) must go over to equation (1). Hence 
A,=7'/l, and y, is the function which represents the mode of vibration of the 
ordinary cantilever beam, i.e. the function y in Part I. Also, the equations obtained by 
substituting the second of equations (16) in equation (2) must hold for all values of z 
and £8; hence, each y,, y,, y.,... Satisfies equation (2). 


The work may be simplified by writing =k8, and regarding k as fixed and £ as 
a variable. From (15) we then have 


and from equation (16) A=A,+A,B+A,8? +... 
(18) 
and 
where A,=A,k+A,, A,=A,k?+A,kK+A,, 
Y,=y, 4+), and Y,=y,k*+y,k+y;. 
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Substituting equation (18) in (17) and comparing powers of 8 up to the second 
degree gives 


2 


+ 7B +k)y,”’ . (20) 


2 2 
and) =6Y,”” Y,=:A.y,+A,¥,+ +4 7 3+k) "4 


2 
++ 


6 
(3+ k?) xy,/” - QU 


- 


3 


~ 


Equation (20) is obtained by comparing coefficients of @ and using equation (19). 
Equation (21) follows by comparing coefficients of 8? and using equations (19) and (20). 


Multiply equation (20) by y, and integrate over O<x</. The left hand side 
vanishes since, on integrating by parts and remembering that y, and Y, both satisfy 
equation (2), it is equal to 

l 


0 


For the integrals on the right hand side, use the integrals (1), (4) and (11) of Table II. 


8 4 
Then «ft +5 k). 


An? 2 4 ] 4 
(23) 


To obtain A,, equation (20) must first be solved for Y, subject to the boundary 

conditions that y= Y, satisfies equation (2). By using equation (19), it may be verified 
that the general solution of equation (20) is given by 


A I 


4l 


where A,, B,, C, and D, are constants. The first three terms on the right hand side of 
this expression give a particular integral, and the last four terms supply the comple- 
mentary function. The boundary conditions give 


(0)+C,y,”" O)=9, (0)+B,y,” (0)= 90, 
2+k 


0 
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For these four equations to be consistent, 


A, (0) Yo (/) 


This equation reduces to equation (22) on using equation (14). Also 
5k) =- (1454) (1s 54), 
=- — —— 1 k =~ = 
The following equation may be obtained from equations (21) and (24) 


Y,”” - A, Y,=A.y,+ (34k) A,Ayy, + (174 4k) A xy, + (4+ k)A,x?y, + 


a 


2 
+14, A, + BA, } 7 :) + 


3 
0 


5 


| 2 
+ += k) (9+2k)+ - } x 
2A 2 
+ 


Multiply this equation by y, and integrate over O<x</. The left hand side 
vanishes as with equation (20). The coefficient of D, also vanishes by equation (22 
The remaining integrals are given by the integrals (1), (2), (3), (4), (11), (12), (13), (14). 
(21), (22) and (23) in Table II. Hence, substituting the values for A,, Bs, C.:Wy and A, 
(see equations (22) and (25)) and simplifying the resulting expression, it follows that 


13 13 8 8 
= + 


( 2 a 4 8 8 


The quantities A,, A, and A, are then given, respectively, by the coefficients of k*, k and 
k® in this expression for A,. 
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Buckling of Columns in the Presence of Creep 


SHARAD A. PATEL 
(Polytechnic Institute of Brooklyn) 


yt SUMMARY: This paper is concerned with the solution of the creep buckling of 

(25) columns. Instantaneous elastic and plastic deformations, as well as the transient 

and secondary creep, are considered. Formulae for the critical time at which a 

column fails are presented for integral values of the exponents appearing in 
the creep law. 


1. Introduction 


During the past few years much work has been published on the subject of 
creep buckling. A reasonably complete list of publications can be found in Ref. 3. 
y + In the present work, Hoff’s theory of creep buckling, in which the instantaneous 
elastic and plastic deformations as well as the transient and secondary creep 
” phenomena are taken into account, is used and a general solution for integral 
values of the exponents appearing in the creep law is found. From this solution 
the critical time, that is, the time required for a column to collapse after load 
application, can be calculated. 


NOTATION 


A dot over a symbol denotes differentiation with respect to time. 


A total cross-sectional area of column 
B= 
side C = 


E  Young’s modulus 
A, L length of column 
M _ bending moment 
P axially compressive load 
a amplitude of non-dimensional deflection 
a,, @, particular values of amplitude a 
Acre Critical values of amplitude a 
h_ distance between flanges of idealised H-section 
m, n_ exponents in plastic and viscous terms of creep law 
ind r radius of curvature 


t, ters tim time, critical time and limit time respectively 
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w = 2y/h, non-dimensional deflection parameter 
y lateral deflection of column 
a, = (-1)* [cos (kz/n) 
B, = m!/[(m—-p)!p!] 
strain 

€., € Strains on concave and convex sides of column 
= /(4L’) 
= (-1)""*- [sin (kx /2n) cot (kx /2n) 

A, » coefficients of viscous and plastic terms in creep law 
o stress 

o., 7 stresses on concave and convex sides of column 


o average compressive stress 


2. Creep Law 


The creep law is assumed as 


This relation is the same as that given in Ref. 1. In the present analysis the 
stress « is considered positive when compressive (instead of tensile as in Ref. 1). 
The first term in the right hand member represents the recoverable elastic 
deformations. The second term takes into account the instantaneous non- 
recoverable plastic deformations and also, in an approximate manner, the transient 
creep. The quantities » and m are so chosen that the second term represents the 
intercept of the straight line portion of the secondary creep curve on the strain 
axis less the elastic deformation (see Fig. 1). The symbols m and n stand for 
integers and the values of k, depend on the conditions of loading. 


When o>0, «>0 (loading in compression): 
« >0, «<0 (unloading in compression): k,=0. 
<0, ¢<0 (loading in tension): =(—1)" 
«<0, «>0 (unloading in tension): k,=0. 


Finally, the third term in the right hand member represents the non-recoverable 
creep deformations of the secondary phase. The values of k, are: — 

when « > 0 (compression), k, = 1; 

when o < 0 (tension), k,=(- 1)"*". 
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o/e 


SECONDARY CREEP 


TRANSIENT CREEP | 


Fig. 1. 


Representation of compressive creep curve. 


3. The Column Problem 


If an axially compressive load P is applied to an idealised H-section column, as 
shown in Fig. 2, the stress-load relation can be given as 


2M 
Ah 
2M 
Ah 


(2) 


where o, and o;, respectively, are the stresses in the concave and convex flanges, A 
is the cross-sectional area of the column, / is the distance between the two flanges 
and M is the bending moment at the point under consideration, given by 


in which y is the deflection. Introducing a non-dimensional deflection parameter 
w=2y/h equations (2), together with equation (3), become 


Hence 


(1+ 


o,=o (1 —w) 


(5) 


where =P/A is the average compressive stress. 
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The relation between the strains and the curvature can be given by 


in which <, and <, are respectively strains in the concave and convex flanges of the 
column and r is the radius of curvature. For small deflections, r can be related 


to w by 


1 ‘ 


where x is the axial co-ordinate of a point on the column axis, as shown in Fig. 2, 
and the suffix x indicates differentiation with respect to x. 


Equations (6) and (7) yield 
When the column begins to creep, the compressive strain <, on the concave side 


increases because of the increasing curvature. For this flange of the column the 
deformation law is therefore 


At the same time on the convex side of the column the increasing curvature 


diminishes the absolute value of the compressive strain <,. Hence for this flange 
a; \* 


Equations (4), (5), (9) and (10), together with equation (8), yield 


(- 2( (l+w)"w+ (7) (1) 


Solution of equation (11) will give the deflection-time characteristics for the 
column as long as the strains remain compressive in both flanges. The initial non- 
dimensional deviation w=w,, at time t=0, of the centre line of the column from 
the straight line along which the load P is applied (Fig. 2) can be represented as"':"’ 


w.=a, sin (=). : « (12) 


For the solution of equation (11) the deflection w at any time ¢ is assumed as 


w=asin ; ‘ . (12a) 
128 
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(6) 
the 
ted TENSION FLANGE 
COMPRESSION FLANGE 
(7) 
Fig. 2. 
Column. 
(8) 
ide 
the 
(9) 
Using the collocation method and satisfying the differential equation (11) at the 
- mid-height of the column leads to 
dt _ (1 +.a)"} 03) 
da (1 +a)" -(1-)" : 
”) where <, is the Euler buckling strain 
Integrating, 
n- 
[tan (k=/n) ]” [cot (kz/n)]’~*" (a of "s+ 
2) (a—a,) cot 
+(—1)"*/? tan ltan(k /n)}" x 
1S 
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Xr n 
(* ¢ (m—p)!p! 
Of x 5 —1 for n even, m,’=m-—1 for m even 
=m for m odd 
= et for n odd, m,’=m for m even 


=m-—1 for m odd. 


It follows from equation (4) that the compression changes to tension in the 
middle of the column when a=1. Since equation (13) was derived under the 
assumption of compression in both the flanges, it is valid only for values of a< 1. 
The equation has a singular point at which da/dt increases indefinitely; hence, 
equation (15) is not valid beyond this point. The singularity occurs when the 
numerator of the right hand side of equation (13) vanishes, that is when 


If the value of a=a.,<1 as obtained from equation (16), the lateral velocity of 
the column increases beyond all limits as the deflection a approaches a.,. If the 
value exceeds unity, equation (15) ceases to be valid. In case a., exceeds unity, 
substitution of a, into equation (15) yields an upper limit of the critical time ¢.. ,. 
since compression is assumed to be present in the tension flange throughout the 
buckling process, which reduces the rate of increase in curvature and consequently 
yields a larger time of failure. The time corresponding to a=1 is designated as 
the limit time f;;,,. Upon substitution of a=1 equation (15) reduces to 


Bes Boa 


9 9 = a,") 
a,” + tan? (kz/n) 27 p=1,2,3 P 


0 


(1 —a,) cot 


(p+1)/2 -1 ai 
1 +a, cot? (kz/n) 


my 


p=2,4, 


(— 1)* [cot (kz/n) ] + 


p—2q 


ty, +eot? (kz/n) 
+(- lo |} (17) 


© cot? (kz/n) 
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Beyond this limit, in the tensile part of the convex flange, equation (10) should 
be replaced by 


In the remainder, equation (10) is still valid because the stress there is still 
compressive. Solution of the problem under these conditions is difficult. But if it 
is assumed that equation (18) is valid in the entire tension flange from this limit 
time on, then equation (11) must be replaced by 


Wine =2 + (7) | + 
+ (CY a9 


Again, solution by the collocation method yields 


dt +a)" +(— 1)" (1—a)"J} 


da (1+a)"—(— (1—a)" (20) 


Integration gives the time required for the increase in the amplitude of the non- 
dimensional deflection from a, to a,: — 


(i) For n odd and m even: — 


lo a, ~ a,? + tan? (kt /n) wa 
N k=1,2.3 p [tan (k /”)] q=0,1,2 p—2q ) 
j2 1 +a," cot® (kr / 
+(— DP 08 cot? 


(ii) For n odd and m odd: — 


Cc a,” + tan? B 3 


2B (n—1)/2 m 


1, 


p-1)/2 (_ 
2,8, [tan (kx | = x 


MN k=1,2,3 p= 


3,5 


(a, —a,) cot | 
x (a,?~ 24 —a,?- +(— 1)*? tan (21d) 
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(iii) For n even and m even: — 


x 24 t(— tan~! (a, —a,) tan (kx /2n) (21c) 


1+a.a, tan? (kz/2n)) 


(iv) For n even and m odd: — 


_, (a,—a,)tan(kr/2n) _ 
k=1,3,5 1 +a,a, tan? (kx/2n) 
Bm *3' a,*+cot?(kr/2n) 2B 
~ n k 3,5 log + cot? (kx /2n) n YB [cot (kx /2n)] 
x (kz /2n) (a,”~ 74 —a,?~ 74) + 
q=0,1,2 p-24 
_ Jog tan? (kz/2n) 
+(— 1" 5108 tan? (ke /2n) Gta) 

where [sin (kx /2n) cot (kz /2n). 


Beyond the limiting value of time, for the deflection to reach a,, the substitution 
a,=1 should be made. The singularity of equation (20) corresponds to a displace- 
ment a.;’ which satisfies the relation 


r 


Substitution of a,,’ for a,, with a,=1, in equations (21) yields the time t,4 necessary 
for the lateral deflection amplitude a to increase from unity to a,,’._ Since tension 
is assumed to be present in the entire tension flange throughout the buckling 
process, this adds to the rate of increase of curvature and consequently decreases 
the time of failure. Hence t¢.,,, the lower limit of the time required for buckling, is 


If the initial deflection a, is greater than unity, then tension is present on the 
convex side of the column from the beginning of creep. Equation (20) is valid 
during the entire buckling process in such cases provided it is assumed that the 
entire convex flange is in tension. Equation (15) with the a., value from equation (16) 
gives again the upper limit ¢.,, of the critical time. The lower limit ¢.,, can be 
calculated by using a, for a, and a,,’ for a, in equations (21). 
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Aline Variation of material constants with temperature for 75S-T6 aluminium alloy. 
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ng, is 
4, Example 
(23) 
A column with a slenderness ratio of 99-34 (¢;=0-001) and an initial deflection 
n the a,=0-1 is considered. It is assumed that the column is made of 75S-T6 aluminium 
valid alloy. Further, it is assumed that the applied average compressive stress o is equal 
it the to half the Euler stress. Material constants, E, », m, A and n are found with the 
1 (16) aid of the data presented in Refs. 7 and 8. Fig. 3 shows how these constants vary 
in be with temperature in the range 211-600°F. The results obtained by the use of 
the formulae presented in the paper are given in Table I. 
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TABLE I 
Temperature 11 300 | 375 600 
| 
| 4-02 x 1032 1-93 x 108 101x107 | 2 
2°25 x 10° 4x 104 5x10 


The results clearly indicate that for relatively small variation of temperature 
there may be a significant change in the critical time. For example, when the 
temperature is varied from 300°F. to 375°F. the consequent change in the critical 
time is 22,000 years to 6 weeks. 


5. Conclusion 

The foregoing analysis can be used to find the upper and lower limits of the 
critical time. The critical time is defined as the time at which the magnitude of the 
transverse velocity of the column becomes infinitely large. This differs from the 
analyses of Refs. 4, 5 and 6 where the deflections themselves become infinitely 
large. The infinite velocity occurs at a finite deflection. This analysis may give a 
lower critical time than the non-linear visco-elastic analysis of Ref. 5. If the term 
expressing the plastic deformations is neglected, the solution is reduced to that of 
the non-linear visco-elastic analysis of Ref. 5. It may be noted that in the present 
work it is assumed that the deflection a, at time t=0 is known. If the shape of the 
column before load application is known, a, can be determined in a manner similar 
to that described in Ref. 1. 
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Supersonic Bangs 
Part II 


P. SAMBASIVA RAO 


(Department of Mathematics, The University, Manchester) 


SUMMARY: The non-linear theory of supersonic bangs, obtained in Part I!) for 
a body accelerating along a straight path, is extended to include curved paths. 
The basic theory remains the same. The important parameter, which appears 
in the theory, is the acceleration component along the ray, the rays being lines 
drawn from points on the flight path at an angle cos~! (1/M) with the direction 
of motion. It is found that the only essential effect of the curvature of the 
path is in the modification of this acceleration component to include a term 
due to the transverse acceleration. With this modification the main results 
are formally the same as in Part I. 


The strength of the bow shock is obtained, and it is found that the effect 
of the curvature of the path is more pronounced at points on the inside of the 
curve, and in general it becomes greater as the distance from the body increases. 
A simple asymptotic formula is obtained which predicts the strength of the 
shock with an error of less than five per cent. at distances of the order of a 
hundred body-lengths. Finally, the theory is compared and contrasted with the 

recent work by Warren"). 


1. Introduction 


In Part I of this paper“ the Whitham theory of the shocks produced by a 
body in uniform supersonic flight was extended to include bodies moving with 
variable velocity. However, for simplicity in developing the theory, only the case 
of a straight flight path was considered. In this paper, the necessary modifications 
are obtained for application of the theory to the shocks produced by a body moving 
along a curved path. Also, in this paper, the opportunity is taken to compare this 
theory with the results obtained by Warren“, whose recent paper was published 
after Part I had been written. 


NOTATION 
x position vector of the field point 
X position vector of an element of the body 


X, position vector of the nose of the body 
Received September 1955. 
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distance measured perpendicular to the axis of the body 
distance of an element of the body from the nose 
maximum value of & 

function defining more accurate wavelets 

first zero of F (ny) apart from »=0 

dummy variable 


angular co-ordinate determining the position of a plane round 
the path 


curvature of the path 

time variable 

time when a disturbance leaves the body 
time when the nose of the body crosses the ray 
small increments in 6 and + 

speed of the body 

length of the body 

thickness ratio 

local sound speed 

speed of sound in the undisturbed flow 
Mach number, U/a, 

Mach angle, sin~' (1/M) 

velocity potential 

acceleration component along the ray 
distance along the ray 

expansion ratio 

cross-sectional area of the body 

a,” (M? -1)/A 

defined by equation (3) 


‘ function describing the wave profile 


pressure 
pressure and density in the undisturbed flow 


ratio of specific heats of air 
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2. Basic Theory 


The basic method is unchanged from Part I. The linear (acoustic) theory 
provides a first rough approximation of the geometry of the shocks; the wave fronts, 


ie. the envelopes of sound wavelets, are the linear approximations to the shocks. 


But the linear theory gives only an inaccurate estimate of the shock strengths. In 
fact, the pressure is continuous across a wave front and the estimate of the shock 
strength is taken from the maximum pressure which is attained in the flow behind 
the wave front (Warren). To obtain an accurate evaluation of the shock strengths, 
it is necessary to incorporate certain crucial non-linear features. In the first 
instance, this involves the introduction of the correct propagation speed (equal to 


"the local speed of sound relative to the fluid) for the individual wavelets. Then. 


since the wavelets move with different speeds, they pile up in compression regions to 
form shocks. Thus, shocks are necessarily introduced as an essential part of 
the theory. 


In applying these ideas, it is possible and also desirable to work mainly with the 
geometrical acoustic approximation to the full linear theory; the reasons for this 
were discussed in detail in Part I. In geometrical acoustics an important part is 
played by the rays, which are defined as the orthogonal trajectories of the wave fronts 
and in a sense are the carriers of the disturbance. In this paper the theory is 
developed for a uniform atmosphere, but the method may be extended to other 
cases. It is easily shown that, for a uniform medium, the rays are straight lines 
drawn from each point of the supersonic part of the flight path at an angle 
4n-sin-'(1/M) to the direction of motion, where M is the Mach number of the 


_ body at the point concerned. Thus, the component of the velocity of the body along 
_ a ray is equal to the sound speed. This provides the approximate criterion that an 
| observer subsequently hears a bang if at any time the velocity component towards 
_ him is sonic; for in that case he would be on a ray. The wave front and the rays 
_ are shown in Fig. 1 for the typical case of an accelerating point source. Fig. 1 
| shows that a novel feature of accelerated flow is that the wave front is cusped and 
_ has a rear portion which may be called the rear shock. This is considered later. 
_ For a real body, which is more like a doublet, shocks will arise from both the nose 


and the tail of the body. Thus, if may be expected that the shock system in Fig. 1 
would be duplicated. In practice, however, the rear shock seems to be single, and 


4 the flow pattern is as shown in Fig. 2 (see Ref. 5). 


The theory of geometrical acoustics also provides an approximate determination 
of the amplitude of the disturbance as it moves along a narrow tube formed by 
neighbouring rays. The assumption is made that the energy of the disturbance 
remains constant; that is, energy is not reflected or diffracted. Then, since energy is 
proportional to the square of the amplitude multiplied by the cross-sectional area of 
the ray tube, the amplitude is proportional to {E (s)}~!, where E (s) is the expansion 
ratio at a distance s along the ray tube and is defined as the ratio of the cross- 
sectional area at s to the value at some reference position. More precisely, E (s) 
would be taken as the limiting value of this ratio as the ray tube shrinks up to a 
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single ray. The determination of E is a simple geometrical calculation and the 
details are given in Section 3. It is found that, omitting a constant of proportionality, 
we may take 


2 2 
where a, is the speed of sound in the undisturbed fluid and A is the component of 
acceleration of the body along the ray. Equations (1) and (2) are of the same form 
as for the straight path. The only difference is in the expression for A, which now 
includes a term due to the acceleration perpendicular to the path; thus A depends 
on the curvature of the path and the angular co-ordinate determining the position of 
the plane containing the field point around the path. E obviously tends to zero 
at the cusp (corresponding to s=A in equation (1)), and the approximation of 
geometrical acoustics breaks down. Therefore, near the cusp it is necessary tc 
consider the full linear solution; this is considered in detail in Section 6. For a 
non-uniform medium it would be necessary to include the factor p,a, in the 
expression for the energy, since it would vary with position. 


In geometrical acoustics reflection and diffraction of energy are neglected. But, 
even when this is justified, the theory requires modification to include the non-linear 
effects already mentioned; correct allowance for the dissipation of energy by the 
shocks is essential. It is found that, at large distances, the additional fall-off in 
amplitude due to this effect is represented by a factor {b(s)}~! in the 
amplitude, where 


(3) 


so that at large distances the shock strength varies with s as 


rae 


0 


To obtain this modified theory, a more detailed description of the flow is 
required and it is obtained by taking any flow quantity, the pressure increase for 
example, in the form 


where 7, is the time when the nose of the body is on the ray through the point. The 
function f describes the wave profile and is determined from the initial disturbance 
at the body, while the factor E~! accounts for the change in amplitude as the 
disturbance prupagates along the ray tube. It should be stressed, however, that 
the present unrigorous formulation of this result is used for simplicity and in order 
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to highlight the basic ideas. In fact, it was shown in Part I that, near the wave 
front, the approximate form of the full linear solution is in accord with (5); the 
actual form is 


M> /2 F 


where p, is the density of undisturbed fluid and 
(t-7,- 5). » » 


The F function depends on the shape of the body and is determined in the full linear 
theory. However, if the changes in the direction and magnitude of velocity in the 
time taken by the body to fly its own length are small, the F function is not affected 
by the acceleration and is the same as for uniform motion; it is 


1 


S” (6) dé 


(8) 


where S(£€) is the area of cross section of the body at a distance € from the nose. 
In Section 4 these conclusions are checked briefly for the extension to curved paths. 


For the improved non-linear theory, €, is replaced by », which is now chosen so 


that the wavelets specified by »=constant are determined with the more accurate 


speed of propagation in place of a,. This follows the standard method explained in 


| Ref. 2 and used in Part I. Then, where the wavelets would otherwise overlap, a 
| shock is introduced and its position and strength are determined from the condition 


that the shock speed is the mean of the speeds of the wavelets on either side of the 
shock. The details are set out in Section 5. 


The main result is the expression for the shock strength as a function of s. 
This now depends on the angular co-ordinate round the flight path, the shock being 
stronger on the concave side of the path than on the convex side. Graphs showing 
the results of numerical computations of the shock strength for different curvatures 
of the path and in different planes through the tangent at a point of the path are 
then given. From the graphs it may be noted that, compared with the case of a 
straight flight path, the shock is stronger on the inside of the curve and is weaker on 
the outside. This effect becomes increasingly more pronounced as s increases. Ina 
particular case (Fig. 6), when the curvature of the path is (1 /3000)m.~’, the ratio of 
the shock strength on the inside of the path for 6=4= to that on the outside for 
§=3n varies from 1-1 to 2, as s increases from about 300 to 1,200 body-lengths. 
(The shock strength at the nose is independent of the curvature of the path). For 
larger curvatures this factor can be greater. The effect of the curvature is felt more 
on the inside of the path than on the outside. For example, in Fig. 6, at a distance 
equal to about 1,200 body-lengths, the shock is stronger for 6=47 by a factor 1-6 
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compared with the case of no curvature (?=47); it is weaker for 6=4z only by a 
factor 0:8. As the cusp is approached, of course, the shock on the inside becomes 
much stronger. 


It is shown in Part 1 that, when the acceleration component of the body along a 
ray is positive, the wave front will have a cusp at one point on that ray and there 
will be a rear shock as in Fig. 1. As already noted, geometrical acoustics, and hence 
the correction based thereon, do not apply in the neighbourhood of a cusp. But it 
seems reasonable to suppose that on the rear shock, well beyond the cusp, the theory 
again gives a correct description of the shock decay. It may well be, of course, that 
some energy is lost in the neighbourhood of a cusp, and that in continuing across 
the cusp the shock strength may be appreciably reduced. Nevertheless, it is felt 
that at least the variation of strength of the shock with s will be given correctly. 
Under these assumptions the pressure rise at the rear shock can be estimated, and 
the details are given in Section 6. An asymptotic formula for the shock strength 
for large distances is deduced, and it is found that as s tends to infinity the shock 
strength decays like s~' (log s)-?, which is typical of spherical shocks.“ 


Finally, Section 6 is concluded with a study of the neighbourhood of a cusp, 
working with the full linear theory. In the full linear theory the disturbance does 
not become infinite at a cusp; it is only the approximation of geometrical acoustics 
which breaks down. By incorporating some of the non-linear effects a rough 
estimate of the pressure to be expected near a cusp is suggested. Of course, it is 
not clear what the detailed motion in the neighbourhood of a cusp is really like. 
For example, from the results given by Lilley and others, the flow pattern appears 
to be as shown in Fig. 2, with the bow shock simply ending in the fluid. In these 
circumstances, evaluation of the shock strength in this region is necessarily of 4 
very tentative nature. 
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Fig. 3. 


3. Expansion Ratios of Ray Tubes 


As explained in Section 2, as the disturbance propagates along a ray, its 
amplitude varies like E-+. In this section the expansion ratio E is determined for a 


| general flight path. 


In Fig. 3, A and B are the positions of the body on the flight path at times - 


| and 7+4r respectively. Consider a plane through AB inclined at an angle @ to the 


instantaneous plane of motion at A. In this plane the rays at A and B are 
AC’ and BB’, and AN and BN are the intersections of this plane with the normal 
planes to the path at A and B. If « is the curvature of the path at A then 
AN cos 6=(1/x). The arc AB=Usé;, where U is the speed of the body; hence the 
angle ANB=U5r x cos 6. If u is the Mach angle sin~'(a,/U), the rays AC’ and BB’ 
make angles 4x - » and 4x — » — 6u with tangents at A and B respectively. Therefore 
the angle BOA between the rays is Usz « cos 6- 5u, and at distance s= AC’ along the 
ray, the distance between the rays is 


where BC is the initial distance. 


Since BC=AB cos cos 1, (10) 
it follows that BC’ = [u (ux cos 6 — | br. 
dr 


At a distance s, the cross-sectional area of the ray tube, bounded by the cones of 
rays at A and B and the planes making angles 6 and +44 with the instantaneous 
plane of motion, is given by 56 B’C’s cos ». From (11) this is 


SOS [ U cos (Ux cos | 56 Br, (12) 
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Fig. 4. 


Omitting a factor independent of s, we may now take 


AY 
where \= [ 
COS U cos d- 


Writing M = =, A =a," (M*— cos 6+ 


2 
relation (14) gives A= % 


The components of acceleration alon 
a,dM | dr and @,*M?x; whence it is ea 


(13) 


(14) 


(15) 


(16) 


g and perpendicular to the path are respectively 


sily seen that A defined in (15) is the component 
of acceleration along the ray. 


Finally, the amplitude is Proportional to 
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- Comparing (17) with the corresponding result in Part I, it is seen that they are of the 

' same form and that the curvature of the path only modifies the expression for the 

acceleration component A. In the plane for which 6=47 the curvature has no 

' effect; and the amplitude has its maximum value for 6=0 and its minimum value 

for @==. It is observed from (17) that, when A > 0, and consequently A > 0, the 

| amplitude tends to infinity as s tends to A; this corresponds to the cusp of the wave 

| front on this ray. Since A varies with 6, the occurrence of a cusp on any ray 

depends on the value of 6. In particular, if the speed of the body is constant, a 

cusp occurs only on the rays directed towards the inside of the curve, that is, those 

rays for which |@|< 4x. Fig. 4 shows the wave front, with the rays, for a source 

which moves with Mach number 2 from a straight path into a curved one. The 

_ cusp is seen to occur on the rays drawn towards the inside of the path, while there 
_ isnone on the rays drawn outwards. 


4, Linear Theory 


Following the procedure outlined in Section 2, a check is now made that, in the 
neighbourhood of the wave front, the flow quantities can be taken proportional to 
EF (€,), with €, =U (t-7,-s/a,) in the linear theory. Furthermore, it was found 
in Part I that F(€,) is given by 


& 


[ S” () dé 


0 


where § (€) is the area of cross section of the body at a distance € from the nose; in 

_ this section it is also verified that this remains true even for the general flight path. 

(13) [| The fact that the amplitude agrees with the result of geometrical acoustics shows that 
| the assumptions of geometrical acoustics are valid near the wave front. 


This paper is concerned only with the description of the flow at reasonably 
large distances from the body. The flow is represented by a source distribution 
(15) — along the path of the body. This will be sufficient for slender bodies, but for an 

| aircraft of arbitrary shape a doublet distribution may also be required. 


(16) Let the nose of the body be at X, at time 7 and let U be its speed. Then an 
; element of the path | dX | at the point X is occupied at time 7 by an element of the 
vely body at a distance from its nose, where 
ent 
€=|X-X, (7) |, ‘ (18) 
17) and, since at distance £€ the cross-sectional area is expanding at the rate S’ (¢), the 


| appropriate source strength at the point X is US’ (£) | dX |. 
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Hence at the point x at time / the velocity potential of the flow is given by 


(19) 
the integration being over all those points X on the path from which the disturbance 
can reach the point x at time ¢. Accordingly = is given by 


Assuming that the changes in magnitude and direction of the velocity of the body in 
the time taken to fly its own length are small, the velocity and the distance from the 
point x of any element of the body in (19) can be replaced by the corresponding 
values for the nose when the nose was on the ray through the point x. If the nose 
was on the ray at time 7,, then (19) becomes 


U 
4n | x — X, (7,) | 


To obtain » in the form given in (5), the variable of integration is changed to £ 
and this can be done conveniently by using = as an intermediate variable. Now 
dX/dz differs from the velocity of the body only due to the change of € with -. 
However, it may be shown as in Part I that d&/dr is relatively small and we 
have approximately 


}dX|=UG)d 


The variables © and = are related by (18) and (20). It was pointed out in Part | 
that the major part of the contribution to the integral in (21) is made by the super- 
sonic part of the flight path. Among the elements of the body that contribute to the 
disturbance at the point x at time ¢, the one farthest from the nose, i.e. with the 
greatest ¢. is at the foot of the ray through x. Let €, be this maximum value of £. 
If the point x is chosen on the wave front at time ¢, the corresponding value of £, is 
zero; hence €, is a measure of the distance behind the wave front. Let -, be the 
time at which the element, at distance ¢, from the nose, is on the ray through the 
pointx. Expanding ¢ as Taylor series near 7 =7,, we obtain 


(23) 


The higher order terms are neglected, as in Part I, under the assumption, which is 
true near the shock, that €,/|x-X, | is small. Differentiating (23) it is seen that 


dé 
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by From (21), (22) and (24), the expression for ¢ becomes 


(19) _U? 

urbance 
where s=|s|=|x—X, (r,) |. 26) 


The relation between € and - is given from (18) and (20) as 


(20 
onding where i is the unit vector along the flight path; from this, — ¢, is found to be 
1e Nose ‘ 
(M* —-1)—MA 
(28) 
Q1) where A is the acceleration component defined in (15). (In deriving (28), use is 
7 again made of the assumption that €,/s and the change of direction in a distance €, 
are both small). Substituting from (28) for €,, the potential function ¢ is finally 
d to & written in the form 
Now 
{2«2-1s(1- *)} 
(22) | where AX =a,” (M? -1)/A, as before, and M is the Mach number at time -,. It 


remains to determine €, in terms of ft, 7,, and s. Now, €, is the distance travelled by 


Part | the body during the time (7,-7,); it is therefore approximately equal to 
super. (7, 7») and, since (t- 7,)a,=s, we have 
to the 
h the = 
of £,~U (1 5). . G0) 
| pe From (29) and (30) it is found, neglecting terms of order €,/s, that 
h the si é 
{2(m2—1)s(1- 
(E,) 
ch is {20 Ns (1- 
that 
1 
a4) where F(é,)= (33) 
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Thus the amplitude is given in accordance with geometrical acoustics and the F 
function is the same as in Part I. 


Perhaps it is worthwhile summarising the co-ordinates now adopted. The 
position x of a point is specified by the “ ray co-ordinates ” -,, 6 and s; +, determines 
the point on the path through which the ray is drawn (this point being the position 
of the nose at time 7,), and 4 determines the plane around the path in which 
the ray lies, while s gives the distance along the ray. In terms of these co-ordinates 
the particle velocity and pressure at time f¢ are given by (30), (31) and (32). The 
expression in (31) is the component of the velocity of the fluid along the ray. In the 
region near the wave front, the component tangential to the wave front is small 
compared to ¢,, since across the wave front the tangential component is continuous; 
hence, in this region the flow is approximately along the rays. Thus the flow along 
each ray tube can be considered independently, and in each ray tube the flow is 
similar to propagation in a duct of variable cross section E (s). 


5. Non-Linear Theory 


The expressions in (31) and (32) for the particle velocity 9, and the excess 
pressure p — p, only differ from the corresponding formulae of Part I in the definition 
of A. Now, from (15) and (16), A is given by 


A= 


a,” (M? —1) 

a, dM)" 

2 
[ 1)! Mx cos 6+ 
thus, a term due to the acceleration perpendicular to the path is included. Thus the 
non-linear improvement of these results is exactly as in Part I, and it is unnecessary 
to reproduce the details here. Introducing the excess pressure ratio (p- p,)/p,. 
since yp,=/P,@,", where y is the ratio of specific heats, the flow quantities are 
given by 


Po (1- 
A 
where the wavelets, defined by »(t,s)=constant, are now determined from the 
accurate speed of propagation a+¢,, a being the local speed of sound. Only the 


dependence of » on f and s is shown explicitly, but it also depends on the parameters* 
7 and @ defining the ray. 


*The suffix of +, is dropped hereafter. 
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Fig. 5. 
The relation determining 7 is found to be 
U (t-—7r-s/a,)+ (8 (M?-1)}? F b(s)=4 
where b(s)= —, . (6) 


Evaluating the integral for b, we have 


A>0 


G7) 
A<O. 
Finally, it is necessary to know the relation between s and » at the shock, i.e. the 
distance along the ray at which the wavelet »=constant runs into the shock. For, 
when this is known, equation (35) gives the position of the shock at time ¢ and (34) 
gives the pressure behind it. If the shock is moving into the undisturbed fluid, the 
required relation between s and the value of » at the shock is (see Part I) 


dn 


(38) 
{32 (M?— 1) } F? 


Figure 5 shows the F-curve of an axi-symmetric body of length /, thickness ratio 
8 and shape given by 


{1- (1-7) }. <1, 69) 


where x is measured from the nose along the axis and r is measured perpendicular to 
the axis. It is clear from (38) that large values of s (when b is large) correspond to 
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values of » near to »,, the zero of F (y) (when there is more than one zero, 7, is the 
first one apart from »=0). This leads to a simple asymptotic formula valid for 
large s. For, from (38), 


"a 


| 


as 25/4 (M2 — 1)! (40) 
F (n) = (y+ 1) M7/4 b(s) | 
hence the shock strength for large s is given by 
9 
| F dn’ 
P— Po ‘M3! 0 (41) 


It should be noted that in this formula only a single numerical factor depends on 
the shape of the body. For the body defined in (39) with 6=0-06 (as in Part 1) 


"0 


| F dy’ 


0 


Thus, accepting this as a typical value and taking y= 1-4, (41) becomes 
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TABLE | 
| (iii) iv) (vi) 
"Straight path urved path with constant speed 
Curve Value of Radius of Radius of 
No Ain Acceleration curvature curvature Radius of 
| body-lengths in sec." 3,000 m2. | 1,500 curvature 
Value of@ | Value of inm. 
136 0° | 60° 3,000. 
(2) 1829°6 68 60° 75° 31’ 6,000 
(3) 3659°2 3:4 82° 49’ 12,000 
(4) 0 | 90° 90° 
(5) —3659°2 | — 34 104° 29’ 97° 11’ — 12,000 
(6) — 1829-6 — 68 120° | 104° 29’ — 6,000 
(7) — 9148 | — 13-6 180° 120° — 3,000 


When s is much smaller than A (although still large enough for the asymptotic 
formula to apply), (36) and (41) show that the shock strength falls off like s~*’*; this 
result is typical of cylindrical waves and it was also found in the case of zero 
acceleration. But, in contrast, when from (36) b(s) ~ A! log (s) 
(see Section 6) and we have the result 


P= Pog (logs), 


which is typical of the decay of spherical shocks (see Ref. 6); this final decay is a 
special feature of accelerated motion. 


It is of interest now to compare formula (42) with the one suggested by 
Warren and to consider his results from the present point of view. Warren also 
starts with the linear theory and he uses a formula which is essentially (34) with a 
typical value inserted for F (y). Then, for the special case of uniform velocity, he 
compares this result with the accurate value obtained by Whitham and finds that a 
further factor (M? - 1)'‘M-*/*s~! is required to take account of the non-linear effects. 
He then retains this factor in the formula for the accelerating body. It can be seen 
from (42) that the correct factor is (M?-1)'M~-"/4[b(s)]-'. Thus the effect of 
acceleration appears not only in the original linear formula (34) but also through 
b(s) in the modifying factor; Warren’s formula does not have this latter effect. For 
small s/| |, it is true that b(s) is proportional to s!, so that Warren’s formula is 
initially correct. Ultimately, however, b (s) is proportional to log s, leading to the 
law of decay of spherical shocks, whereas Warren’s formula leads to 


for large s. ahs 
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Calculations of the variation of shock strength with s have been made for the 
body defined in (39) in the case when M=1-05. The strength has been computed 
for various values of A, both from the exact result given by (34) and (38) and from 
the asymptotic formula (42). The results of the exact calculations are shown in 
Fig. 6; the values of s are measured as multiples of body-length / and the values 
of A (again measured in body-lengths) are shown in column (ii) of Table I. 


In Fig. 7 the values of shock strength obtained from the asymptotic and actual 
formulae are compared for the values of A equal to 1829°6, and - 9148, 
measured in body-lengths. It should be noted that the error, which is about five 
per cent., does not vary much over the range of s shown. This is because » 

q 
approaches 1, only slowly and consequently | F (') dy’ is approximately constant 
over an appreciable range. although it is still not equal to { F (1') diy’. 


0 


Thus, even at distances of the order of a hundred body-lengths, the law of 
variation with s of the shock strength is given quite accurately by the asymptotic 
formula, but the “ constant of proportionality ” is to be reduced. 


Now. A is determined by the component of acceleration of the body along the 
ray, and this may be derived from various combinations of the accelerations along 
and normal to the path. Thus, for example, a given value of A may be interpreted 
as corresponding to a body accelerating along a straight path or to a body moving 
with constant speed on a curved path. 


In columns (iii) to (vi) of Table I, various possible interpretations of the set of 
values of A have been given. 
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Since the sound speed a, enters into the expression for A, it is necessary for 


| these results to specify a definite value for the body-length; here / is chosen to be 
_ one metre. In column (iii) the values of A are interpreted as different accelerations 
| of the body along a straight path. In columns (iv) and (v) they are associated with 
| different angles @ about the flight path for a body moving with constant speed on a 
4 curved path; in column (iv) the path has a radius of curvature of 3,000 m. while in 
- column (v) it is 1,500 m. The infinite value of A corresponds to uniform velocity 


and also to the plane 6=4x. The rays for which | |< 4x correspond to the 
inside of the path and the rays for which |4|> 4x correspond to the 
outside. It may be noted from Fig. 6 that the effect of the curvature is more 
pronounced on the inside than on the outside of the path. For example, when the 
radius of curvature of the path is 3,000 m., at a distance of 1,200 m. from the nose, 
the shock is stronger for 6=4= than for @=4= by a factor 1-6, while it is weaker for 
@=2z only by a factor 0-8. Finally, in column (vi) the values of A are interpreted 
for points in the plane of flight for different values of the curvature. 


_ 6. Rear Shock and Neighbourhood of a Cusp 


This section deals with the application of the theory to the rear shock and 
indicates a method of obtaining an estimate of the shock strength near a cusp. It 
has already been mentioned that. if the acceleration component along a ray is 
positive, the wave front has a cusp on the ray and there is a rear shock. On the 
rear shock, well beyond the cusp, it is expected that the geometrical acoustics, as 
modified in Section 5, can again be applied. It is therefore assumed that the shock 
strength decays with s according to the asymptotic formula (41). The formula for 
E (s) given in (13) becomes negative beyond the cusp (s > A) and, of course, the sign 
must be changed; similarly for b (s) we now have 


0 


when s >A. Hence, the asymptotic formula becomes 


"0 


4 
| F di 
(y+1)! A! [47 + cosh~ '(s/A)!] s(s — A) 


However, it may well be that the “ constant of proportionality ” will now be smaller. 
since there may be a loss of energy at the cusp. Certainly. experiments seem to 
indicate that the rear shock is always very much weaker than the main shocks. 


For very large s (s > >A), we have 


{| F di’ } 


P- Py _ 23/4 Me 'ya,! 
(y+ DEA! s (log s) 4 
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This form applies very soon at the point where the rear shock meets the axis, 
since A tends to zero there. Again, for very large s, formula (47) indicates a decay 
of shock strength typical of spherical shocks. 


Figure 8 shows the shock strengths plotted against s for the front and rear 
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' shocks for A= 914-8, 1829-6, and 3659-2, the unit of length being the body length. 
» The shock strength is seen to fall off very rapidly on the rear shock, so that this 
_ again predicts that the rear shock will be very much weaker than the main shock. 


As mentioned in Section 1. geometrical acoustics, and hence the corrected 
theory, break down near a cusp. But the full linear theory gives a finite value for 
the pressure rise and, if a method can be found of modifying this to include the 
non-linear effects, some estimate of the shock strength near a cusp can be obtained. 


Figure 9 is a typical graph of the pressure rise plotted against 
£,=U(t-7r-s/a,), (ie. U times the time that has elapsed since the wave front 
_ passed the point). Now, in the linear theory £,=0 corresponds to the wave front 
and the pressure is continuous there. But in the improved theory, which is applied 
_ away from the cusp, the appropriate value of €, is simply the 1 determined from (38). 
' The portion of the pressure curve (0<£€,<y), which is cut out in this way, 
» corresponds to the wavelets which have been absorbed by the shock. Near the 
' cusp, where the approximate formula of geometrical acoustics breaks down, the 
| appropriate pressure curve must be obtained from the full linear theory. (The 
| pressure is again plotted against €,—U(t-7-s/a,)). Then, having found the 
pressure curve, the problem reduces to choosing the appropriate point on the curve 
_ which corresponds to the value just behind the shock. 


Now, since wavelets are absorbed by the shock at a relatively slow rate, it seems 
| reasonable to neglect this effect just in the neighbourhood of the cusp, and to 
' determine the required point on the pressure curve from the value of €, given by the 
_ theory away from the cusp. In fact, since the integral defining b (s) is convergent as 
s tends to A, the value of » determined by (38) should still be taken. In many cases 
i considered the cusp will be at a large distance, and the approximate determination 
, with » near to n, may then be used. Since, for points in the neighbourhood of a 
cusp, the exact linear expression for p- p, does not vanish for »=»,, it is perhaps 
' possible simply to choose 7 equal to , for much of this region. 
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The expression for ¢ in the full linear theory has been given in Section 4 but 
the formula for the pressure increase p-p,=-po? is now required. The 
disturbance is represented by a distribution of sources on the flight path and at any 
point X the source strength is US’ (€)| dX |, where € measures the distance of X 
from the nose (i.e. the point X, (7) ) at time z. The contribution of this source to the 
value of ¢ at the point x at time ¢ is 


US’ (§) | dX | 


4x|x-X| ’ 
where the retarded time r, used in the determination of €, satisfies 


Clearly, the rate of change of this contribution with time is due to the rate of change 
of S’ (€) as the body passes the point X and this is US” (€). Hence, the contribution 
of the elementary source to ¢, is 


and ¢, is given by the integral of this quantity over all the sources. Making similar 
approximations used in (21) gives 


Then, using the approximation | dX | ~ U dr (see (22) ), we write 


” (2) de 
Po 


where =| X - X, (7) | and X is determined in terms of 7 by (48). 


The method just suggested. although possible in principle, would be cumbersome 
in practice unless some approximation going beyond geometrical acoustics can be 
used for (50). This sort of difficulty arises in many applications of geometrical 
optics and acoustics; near certain curves, called “caustics,” the theory of 
geometrical acoustics is invalid and a more accurate approximation is required. 
Methods for overcoming these difficulties have been developed for some of these 
cases and it may be that similar ones may be applied here. On the other hand. it 
may be possible to obtain adequate results by retaining further terms in (23) (the 
second order term vanishes at a cusp). However, it is felt that such investigations 
would be more appropriate when a reasonable qualitative account of what happens 
in the neighbourhood of cusps is available. 


In Warren’s® approach to the problem, it is not necessary to have such an 
accurate estimate of (50) and he obtains a rough formula as follows. He suggests 
a method for estimating the shock strength at the cusp which is derived from a 
formula corresponding to (50). From the present point of view it may be explained 
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in the following way. First, the integral (50) is approximated by some mean value 
of S” (€) multiplied by Ar, where Az is the change in the retarded time 7 over the 
range of integration. At a fixed point x, Ar varies with the time ¢ and Warren 
chooses the maximum value, giving the maximum pressure rise as the linear estimate 
of the shock strength. Then, as before, he introduces the factor (M? — 1)'M~7/*s~* to 
account for the extra decrease in strength due to non-linear effects. In contrast, in 
the method explained, this extra decrease from the linear result arises because, not 
only is the whole pressure curve of Fig. 9 diminished as s is increased, but the point 
on the pressure curve which represents conditions behind the shock tends to the 
zero Of p- p, as s tends to infinity. 
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The expression for ¢ in the full linear theory has been given in Section 4 buf ; 
the formula for the pressure increase p-p,=-p,.%: is now required. Th ® 
disturbance is represented by a distribution of sources on the flight path and at any} 
point X the source strength is US’ (€) | dX |, where € measures the distance of Xf 
from the nose (i.e. the point X, (r)) at time r. The contribution of this source to the 4 


value of ¢ at the point x at time f is 


US’ (&) | dX | 


’ 
where the retarded time r, used in the determination of £, satisfies 


Clearly, the rate of change of this contribution with time is due to the rate of change 
of S’ (£) as the body passes the point X and this is US” (é). Hence, the contribution 
of the elementary source to 4; is 


| dX | 


and ¢, is given by the integral of this quantity over all the sources. Making similar 
approximations used in (21) gives 
P- Po _ _ yM? 


Then, using the approximation | dX | ~ U dr (see (22) ), we write 


” () dr 5 


where =| X - X, (r) | and X is determined in terms of = by (48). 


The method just suggested, although possible in principle, would be cumbersome 
in practice unless some approximation going beyond geometrical acoustics can be 
used for (50). This sort of difficulty arises in many applications of geometrical 
optics and acoustics; near certain curves, called “caustics,” the theory of 
geometrical acoustics is invalid and a more accurate approximation is required. 
Methods for overcoming these difficulties have been developed for some of these 
cases and it may be that similar ones may be applied here. On the other hand. it 
may be possible to obtain adequate results by retaining further terms in (23) (the 


second order term vanishes at a cusp). However, it is felt that such investigations | 


would be more appropriate when a reasonable qualitative account of what happens 
in the neighbourhood of cusps is available. 


In Warren’s® approach to the problem, it is not necessary to have such an 
accurate estimate of (50) and he obtains a rough formula as follows. He suggests 
a method for estimating the shock strength at the cusp which is derived from a 
formula corresponding to (50). From the present point of view it may be explained 
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' in the following way. First, the integral (50) is approximated by some mean value 
of S” (€) multiplied by Ar, where Az is the change in the retarded time = over the 
' range of integration. At a fixed point x, Ar varies with the time ¢ and Warren 


chooses the maximum value, giving the maximum pressure rise as the linear estimate 
of the shock strength. Then, as before, he introduces the factor (M? — 1)'M~*/*s~! to 
account for the extra decrease in strength due to non-linear effects. In contrast, in 
the method explained, this extra decrease from the linear result arises because, not 
only is the whole pressure curve of Fig. 9 diminished as s is increased, but the point 
on the pressure curve which represents conditions behind the shock tends to the 
zero of p— p, as s tends to infinity. 
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The Behaviour of Damped Linear Systems in 
Steady Oscillation 


R. E. D. BISHOP 


(University Engineering Laboratory, Cambridge) 


Summary: The classical theory of small harmonic vibrations of a linear 
damped system embodies the notion of “ viscous damping.” The equations of 
motion which result are somewhat complicated and, when there are more than 
two degrees of freedom, they are usually too unwieldy to be of much practical 
value. When the damping is small, however, approximating assumptions 
may be made which permit the treatment of systems which are near 
resonance as if they possess but one degree of freedom."!) But the effects 
of making these assumptions are by no means easily assessed, and even 
their justification is tedious. 


It is shown that these difficulties may be greatly diminished by postu- 
lating hysteretic damping instead of viscous damping; the concept of 
hysteretic damping has been dealt with in two previous papers.) “) The 
equations then take a much simpler form and the justification for, and validity 
of, the foregoing approximations are more easily seen. Moreover, the 
effects of damping upon the principal modes which the system possesses in 

the absence of its damping may be elucidated in this way. 


1. Introduction 


When a mechanical system having a finite number of degrees of freedom 
executes steady forced harmonic vibration, the motion that the system would 
perform if it were conservative is modified by the existence of damping forces. If 
the excitation frequency is remote from a natural frequency of the undamped 
system, this modification is not great unless the damping is very heavy; but the 
damping forces play a dominant part in determining the amplitude of the response 
when the impressed frequency is approximately equal to a natural frequency. The 
engineer is most concerned with this latter problem, particularly when the damping 
is light, since it is intimately connected with the liability of components and 
structures to fracture by fatigue and with the other objectionable features of violent 
oscillation. 


The usual way of tackling the problem of small dissipation is to postulate the 
existence of viscous damping forces. This leads to a theory which is very 
complicated when applied to any but the simplest systems. As a result, practical 
vibration analysis”? usually involves the assumption that the system may be 
thought of as possessing one degree of freedom when it oscillates at, or near. 
resonance. That is to say, it is supposed that the response at the principal 
co-ordinates of the system (as they are defined in the absence of damping) may be 
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analysed independently such that, at each resonance, one of these co-ordinates only 
varies. The maximum response in each of the modes is assumed to occur when the 
excitation frequency equals a natural frequency of the system (which is calculated in 
the absence of its damping) and the variations of the amplitudes of the principal 
co-ordinates with frequency is supposed to be the same as that of the displacement 
of a simple damped spring-mass system. 


This method of approach is justified when the damping is light. But in fact 
there is always coupling between the modes, although one of them predominates at 
each resonance. Also the maximum responses do occur at frequencies which are 
very Close to the natural frequencies. These results are not easily demonstrated 
with any generality and an algebraic discussion of their consequences in a given 
system is usually out of the question as a result of mathematical complexity. 


The purpose of this paper is to show that the qualitative results just mentioned 
may be arrived at with comparative ease by changing the damping theory to that of 
hysteretic damping®:*. Moreover, the algebra is materially simplified by this 
change of hypothesis and provides greater insight into the nature of the mathe- 
matical theory of damped vibration; it also greatly simplifies the analysis of systems 
whose damping is not small. The method which is used to demonstrate this is to 
show how the introduction of damping modifies certain of the simple results of the 
standard theory for conservative systems. The first part of this paper is therefore 
devoted to a brief summary of the relevant features of the theory of 
undamped systems. 


The reader should note that it is not the author’s purpose to present a practical 
technique of vibration analysis for use with known systems. It is, rather, to discuss 
the nature of damped vibration in general terms. 


NOTATION 


real constants defined in equation (14) 
complex constants defined in equation (25) 
inertia coefficients for kinetic energy 
stability coefficients for potential energy 
hysteretic damping coefficients 

v¥(-1) 

number of degrees of freedom 
generalised force corresponding to q, 

generalised co-ordinate 

kinetic energy 


potential energy 


receptance function (see equations (7) and (18) ) 
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4 determinant defined in equation (6), but which may have 
Crs replaced by c,,+id;s 


+%n-, complex roots of 4,,=0 
€,,&,....&, complex roots of A=0 
®, amplitude of Q, 
amplitude of q, 
anti-resonance frequencies 
© excitation frequency 


©,,.,...@, natural frequencies 


2. The Conservative System 


2.1. THE GENERAL SOLUTION 


The kinetic and potential energies of a linear oscillatory system with n degrees 
of freedom have the form 


T =4 [@,,9,7 + AnnGn® +. ] (1) 


V =4 + +. Can Gn" + ] 


If these functions are used in conjunction with Lagrange’s equations, the general 
equations of motion are found; namely, 


where r=1,2,. . .n. Suppose that the excitation is harmonic such that 
O,= ®, el" n) . (3) 


and that the harmonic response is sought in the form 


Equation (2) now gives, for the amplitudes Y,, 


(c,,-— ©74,,) + (Cop + (Car — V, (5) 
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To solve this set of n equations, it is convenient to form the determinant 


2 2 2 
H 9 2 

| 
| | 
Cay An, Cno Ans . Con— Ann | 


Suppose first that Q, is the only non-zero generalised disturbing force; the 
amplitude of g, is then given by 


where A,, is that determinant which is formed from A by omitting the r‘" column 
and s““ row. In general all the ®’s will be non-zero and, in that case, the Y’s must 
be found by adding solutions of the type (7). 


The quantities z,, have been called “ receptances””; they are functions of the 
energy coefficients of equations (1) and of ? and they have been shown by Duncan“? 
and others®: © to possess properties which are both interesting and useful. It may 
be mentioned that these same properties are still available when linear damping is 
incorporated in the theory. Attention will now be paid to the algebraic form of 
these receptances. 


The natural frequencies of the system may be regarded as those excitation 
frequencies for which an infinite response is obtained from a finite amplitude of 
disturbance. They are thus given by the roots ,’, ,”,..... .,” (in order of 
increasing magnitude) of the frequency equation 


The quantity A is a polynomial of degree n in ? and the roots of (8) are all real 
and positive provided that the system is a stable one; mathematically, this is because 
the expressions (1) for T and V are “ positive definite.’ The receptance z,, may 
therefore be written in the form 


__ fe’ 


| where f (w?) is a polynomial of maximum degree n-1 in w*. This general result 


' will now be used to demonstrate two general results of vibration theory. 
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2.2. DIRECT RECEPTANCES 


Suppose that it is required to measure the response at qg, caused by excitation 
Q,; the receptance of equation (7) would be denoted by z,, and is referred to as “ the 
direct receptance” at q,. It is given by 


(10) 


Now the equation that is obtained by setting A,, equal to zero is the frequency 

equation of that system which is obtained from the original one by locking it in such f 
a way as only to prevent variation of g, without adding mass or elasticity at any 
point. It follows that the equation 


has n-1 roots which are all real and positive and they will be denoted by 
0,7, 2,7,.... O2,_,7.. Thus a direct receptance may be written, not only in the F 
form (9) but more particularly as 


(12) 


where K is a constant. The frequencies , are the natural frequencies and the °2, 
are the so-called “ anti-resonance frequencies.” When © coincides with one of Ff 
these anti-resonance frequencies, no motion occurs at the driving point q,. 


The anti-resonance frequencies of a mechanical system are spaced between the 
natural frequencies in order of magnitude. With this limitation, values of these 
two sets of frequencies may vary between limits which are set by the possibility of 
having repeated roots of equations (8) and (11). These special limiting cases (which 
of course, cannot be realised exactly in a practical system) have been examined by 
Bishop and Johnson. 


2.3. MOTION IN PRINCIPAL MODES 


A second property of the receptances 2z,, of equation (9) is that all of them may 
be expressed as series of partial fractions. Thus z,, may be expressed in the form 


(13) 


where the numerators are all real constants. 


The significance of this expansion can easily be seen by writing down the 
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responses at the co-ordinates q,, q2, .-- qu due to harmonic forcing at one 
co-ordinate only—the s‘" say. They are 


— 


— w? @? 


If excitation is applied at more co-ordinates than just the s‘", then other expressions 
of the same form as these must be added, but the general form of the solutions 
is retained. 


Each column of terms in equation (14) represents motion in a principal mode. 
Let » approach the m' natural frequency ,, and, at the same time, “, be 
diminished such that the fraction ,/(«,,” — *) retains a finite and constant value. 
As this limiting process proceeds, the harmonic distortion acquires a definite shape 
which is given by 


This shape of distortion is, of course, independent of the co-ordinate g, at which the 
vanishingly small excitation ?,e*” is applied; in fact, the appearance of the suffix s 
in the relation (15) merely indicates that, while the ratios between the A’s are 
determined by the system, their actual values depend upon the co-ordinate gq, at 
which forcing occurs. 


When w does not coincide with a natural frequency, equation (14) shows that, in 
general, all the principal modes are excited. Moreover, if a principal co-ordinate p,, 
is associated with distortion in the m™ principal mode—that is with distortion 
according to the m' column of partial fractions—then p,, varies in the same way 
with w as does the deflection of a simple spring-suspended mass when it is subjected 
to a harmonic force. 


2.4. THE ADDITION OF DAMPING 


The various results that have been quoted for conservative systems are well 
known. They are modified when damping is included and it is difficult to see the 
ways in which they are altered by viscous damping, because of the algebraic 
complications which accompany viscous damping terms. Now heavier algebra is 
evidently inescapable when damping is introduced, but it will next be shown that the 
theory does not lose its previous form and indeed retains much of its simplicity 
when the existence of hysteretic (instead of viscous) damping is postulated. 
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3. Hysteretic Damping of Systems with m Degrees of Freedom 
3.1. THE GENERAL SOLUTION 
The mathematical theory of hysteretic damping has been presented in a 


previous paper“, and it will now be supposed that all the damping forces are of this 
type. It was shown that the equations of motion governing the variation of 


+ + (Cor +ider) (16) 
t+ (Cue + G.] =P, (p=1,2,...a) 


where the effects of the hysteretic damping are embodied entirely in the constants 
d,,,d,.. . .,and where the applied harmonic forces give rise to generalised forces, 
as in equation (3). The algebraic consequence of introducing the damping into a 
conservative system is thus to add an imaginary part to each of the stability 
coefficients c,, of the potential energy function (1). If the harmonic solutions of 
equation (4) are now sought as before, n algebraic equations are found for the 
amplitudes V;. They are 


and the same technique of solution as before may be used. 


The amplitude of the response at q, due to excitation by a force P,e is again 
found to be given by equation (7), where the receptance 


is now complex. This means that the response is no longer in phase (or in exact 
anti-phase) with the excitation. The complexity of 2,, arises because the 
determinant A is similar to that of equation (6) but has the real constants c,, of that 
expression augmented by the imaginary parts id,,, and also the determinant 4. 
which is formed from it is complex for the same reason. 


When it is expanded, the determinant A yields a polynomial of degree n in ©” 
whose coefficients are now, in general, complex. Equation (8) will have n complex 
roots, which will be denoted by €,, €.,. . . . &n, So that the receptances have the form 


where f (@) is now a complex polynomial in ? of maximum degree n-1. The 
receptances may be manipulated as before to show (i) the properties of direct 
receptances and (ii) the nature of motion in modes although these modes are no 
longer the “ principal modes.” 
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3.2. DiRECT RECEPTANCES 
The direct receptance 2,, is given by 2,,=A,,/4 and, just as A=0 has n 


' complex roots, so 4,,=0 has n-1 complex roots which will be denoted by 


aS ae Nn-1. The receptance may therefore be written in the form 


_ where K is a constant. 


If all the damping forces are made smaller and smaller then, in this expression 


(21) 


oF; 


in order that the result (12) may be arrived at for the conservative system. 


The quantity «? is essentially real and positive because the motion is restricted 
in this way so that, in general, | z,, | can never be made infinite or zero by adjust- 
ment of the excitation frequency. The amplitude | V, | which is given by 


may be said, somewhat loosely, to have minimum values when 


w=R(n) (i=1,2,...n-1). (23) 


and maximum values when 


where it is assumed that the introduction of damping does not modify the quantities 
Q? and w;? to the extent of making the real parts of 7; and ¢; negative. This 
assumption is justified as normally damped systems do give minimum and maximum 
amplitudes in the manner just stated. 


3.3. MOTION IN “ MODES” 


The receptances z,, have been shown to be of the form (19) where f (w7) is a 
polynomial of maximum degree n-1 in w* whose coefficients are complex. It 
follows that they can still be expressed as series of partial fractions as in 
equation (13) so that 


tn : 


where the constants ,,B,, and (of course) the £’s are complex. 
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As before, it is instructive to write down the expressions giving the responses at 
the co-ordinates g which are caused by harmonic forcing of amount Q,=°,e 


They are 


Sn - 


Once again, if more of the generalised forces than just Q, are non-zero, then other, F 
similar, solutions must be added to these; but they will not alter the general 
character of the equations. However, there is no longer any assurance that the 
ratios between the constants ,,B,, in any column of partial fractions are the same for 
all co-ordinates g, at which excitation is applied, since the previous limiting condi- 
tion cannot now be used. Thus the shapes of the “ modes” which are associated 
with the various columns depend upon the nature of the applied harmonic 
excitation. 


It will be seen from equations (26) that, for any value of , all the modes will 
be excited in general; but one mode, the m™ say, will predominate if 


(27) 


Em) = 0? 
SF (&m)=0 


Also, the amplitude of the motion can never be infinite unless these approximations 
are realised exactly. 


If a particle of mass M is suspended by a light spring of stiffness k which has 
in parallel with it a hysteretic damper whose coefficient is h, the equation governing 
extensional displacement x if a harmonic force Fe! is applied to the particle is 


The steady-state solution of this is 


Feit 
(29) 
(k+ih)- Mo® 
This result has the same form as the contributions of individual partial fractions to 
the variations of the q’s in equations (26). The nature of the approximation which 
is made when a damped system is treated as if it possesses but one degree of 
freedom when it is excited at a near-resonant frequency may thus be seen clearly. 
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Equations (26) provide a means of assessing the errors arising in a 
‘vibration analysis which ignores damping altogether. This latter technique is, of 
‘course, only feasible when the frequency is not near a natural frequency, as will be 

seen on inspection of the equations. When the damping is not light, it will often be 
| possible to use only a few of the partial fractions instead of all of them; the decision 
‘to make a selection of this sort may be made when numerical values of the various 
| parameters are known. 


4. Systems with Small Damping 
'4.1. THE GENERAL SOLUTION 


In engineering one is mainly concerned with systems which embody “ small ” 
| damping, as these are the most prone to violent oscillation near resonances. It is 
| worth while, therefore. to discuss the form that the theory of Section 3 takes when 
' the damping forces are small. It is convenient, in doing this, to quote certain 
; mathematical results, the validity of which is dealt with in the Appendices. 


Suppose that a conservative system is modified by the introduction of small 
_ damping forces, so that the coefficients c;; have, added to them, small imaginary 
parts id;;. The quantities €,, €.,... . &, which are the roots of A=0, may be 
; regarded as functions of the coefficients (c;;+id;;). By means of Taylor’s theorem 
| for a function of several variables, it is shown (see Appendix I) that 


En =,” +i [real function of order (d;;) ] + | (30) 
+ [functions of higher orders of (d;;)]_ (m=1,2, . . . n) 


| The condition under which the first term in square brackets is small and the second 
| is negligible is that of “ small damping,” and its nature is discussed in Appendix I. 
Thus, for small damping, any receptance may be written in the form 


f 


= (w 7+ iO [dis] +iO [d;;} — w*) + iO [dij] (31) 


Once more it is convenient to examine separately (i) the properties of direct 
receptances and (ii) the nature of the motion in modes. 


| 42. Direct RECEPTANCES 
The roots of A,,=0 are modified in a similar fashion to those of A=0, so that 


tn +iO [dij] (m= 1, 2, 1) (32) 


Where the function of order (d,;) is a real one. Direct receptances may therefore be 
Written in the form 


(@,7 + iO [dij] ©?) (,? + iO [dij] ©?) +10 [dij] ©”) ) 
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It is now seen that, for small damping, the frequencies which correspond to the 
minimum motions at g, are very close to the original anti-resonance frequencies of 
the undamped system. Furthermore, these minimum responses are small. Equally, 
maximum motions occur at excitation frequencies which are very close to the natural 
frequencies of the conservative system; these motions are large. 


4.3. MOTION IN MODES 


In Appendix II it is shown that any receptance 2,, of the form (25) is such that, 
for small damping, 


+ iO [dij] (m= 1, 2, n) (34) 


where, again, the function of order (d,;) is real. It follows therefore that. if 
excitation is applied at q, as before, 


iO [dij] 2A as + iO [dis] 


As, + iO [du] 


92 [dy] - 0? 0,7 +i0 [dy] - On? +10 [dij] - 


An +iO Edu], +iO [du] iO 


(35) 

The modes are seen to be still complex but the complexity, that is the amount 
by which the displacements are out of phase in a given mode, is small. The modes 
are very nearly the same as the principal modes but no longer represent motions 
which are exactly in phase. 


5. Conclusion 


The theory of small forced harmonic motions of conservative systems is well- 
known and is mathematically straightforward. But when damping is included in 
the theory, added complication of the equations is inescapable. The usually adopted 
theory is that of viscous damping and, by its use, it is possible to derive several 
important results for the case of small damping. These results are not easily 
established and the general theory of viscous damping is, analytically speaking. 
unwieldy. Problems in which the damping is not small are by no means easy to 
tackle on this account. 


It has been shown in this paper, that the motion of damped systems is more 
easily dealt with under the assumption of hysteretic damping. The same general 
conclusions are arrived at as may be found from viscous damping theory, but their 
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extraction is much simpler. The hysteretic damping theory also provides a simpler 
means of dealing with systems whose damping is not small. 


ACKNOWLEDGMENT 
The author wishes to acknowledge the help which he has received from 
discussions with Professor D. C. Johnson on the subject matter of this paper. 


REFERENCES 


1. Yates, H. G. Vibration Diagnosis in Marine Geared Turbines. Transactions of the 
N.E. Coast Institution of Engineers and Shipbuilders, Vol. 65, 4, 225-261, 1949. 


tv 


BisHop, R. E. D. The Treatment of Damping Forces in Vibration Theory. Journal of 
the Royal Aeronautical Society, November 1955S. 


3. BisHop, R. E. D. The General Theory of “ Hysteretic Damping.” The Aeronautical 
Quarterly, February 1956. 


4. Duncan, W. J. Mechanical Admittances and Their Applications to Oscillation Problems. 
R. & M. 2000, 1947. 


5. JOHNSON, D. C. The Application of Admittance Methods to the Classical Theory of 
Small Oscillations. Engineering, Vol. 171, No. 4453, p. 650, 1951. 


6. BisHop, R. E. D. The Analysis and Synthesis of Vibrating Systems. Journal of the 
Royal Aeronautical Society, October 1954. 


7. BisHop, R. E. D. and JoHNson, D. C. Some Properties of Nodes in Vibrating Systems. 
The Aeronautical Quarterly, February 1955. 


Appendix I 


THe EFFECT OF SMALL DAMPING ON ROOTS OF THE EQUATION A=0 
The equation A=0 . . (36) 
may be regarded as one relating any of its roots w?=€, to the complex quantities 


in which only the imaginary parts will be supposed to vary. Thus, equation (36) may be 
written as 

the quantity €, being given as an implicit function of the z,. The value of €, varies if 


the z,, are varied, and, in particular, if all the imaginary parts of the z,; are zero, 
then £,=w,?. 


According to Taylor’s theorem for several complex variables 


where the derivatives are all evaluated at the real values of the z;;, (that is at £,=w,’). 
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Now, in order to obtain these derivatives, it is necessary to differentiate the implicit 
function (38) of €, and the complex variable z,,. This may be done by means of a 
result in the theory of complex variables which states that, provided that conditions of 
analyticity are met and that equation (36) has no repeated roots, 


For any physical system, the conditions under which this result holds good are fulfilled. 


The two derivatives on the right hand side of equation (40) are real. It follows 
that the one on the left hand side is also real and hence that the second terms of the 
Taylor series (39) are all imaginary, containing the d,; to the first power. 


This result provides a mathematical specification for “ small damping,” namely that 


for all i, j and r. This implies that all the | d,;; | must be much smaller than the absolute 
values | c,, | of these stability coefficients which are non-zero. 


Appendix II 


THE EFFECT OF SMALL DAMPING ON THE CONSTANTS ,B,, 


The constants ,B,, of equation (25) can be expressed in a standard form, as may be 
seen by starting from the equation 


= 2 
A w? (4 ) 
If both sides of this equation are multiplied by A and then w? is made equal to &,, it is 
found that, provided that A=0 has no repeated roots, 


(; 


This equation defines ,B,, as an analytic function of the quantities z,,, (see equation (37) ) 
and the ¢’s. But, since the é’s are analytic functions of the z’s, it follows that the 
;B,, are analytic functions of the z’s. 


Taylor’s theorem may now be used in the same manner as before to trace the 
effect upon the constants ,B,, of varying the damping coefficients d,,,._ It will be recalled 
that, when all the d,, are zero, ,B,, is the same as ;A,, of the conservative system, so that 


u=1 


where all the derivatives are evaluated at one of the values £=w,? (when all the z’s are 
real). From inspection it is clear that the derivatives in the second terms of equation (44) 
are all real so that, to the first power of the d,,, ;A,, is modified by the addition of an 
imaginary term only. 
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